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Wc discuss several pairing-related phenomena in nuclear systems, ranging from superfluidity in 
neutron stars to the gradual breaking of pairs in finite nuclei. We focus on the links between 
many-body pairing as it evolves from the underlying nucleon-nucleon interaction and the eventual 
' experimental and theoretical manifestations of superfluidity in inflnite nuclear matter and of pair- 

' ing in flnite nuclei. We analyse the nature of pair correlations in nuclei and their potential impact 

on nuclear structure experiments. We also describe recent experimental evidence that points to a 
, relation between pairing and phase transitions (or transformations) in finite nuclear systems. Fi- 

nally, we discuss recent investigations of ground-state properties of random two-body interactions 
■ ■ where pairing plays little role although the interactions yield interesting nuclear properties such 

^ ' as 0+ ground states in even-even nuclei. 
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I. INTRODUCTION 

Pairing lies at the heart of nuclear physics and the quantum many-body problem in general. In this review we address 
some of the recent theoretical and experimental studies of pairing phenomena in finite nuclei and nuclear matter. In 
infinitely extended nuclear systems, s uch as neutron star matter and nuclear matter , the study of superfluidity and 



pairing has a long history, se e e.g., (Cooper et a/. , 195£ ; Emery and Sessler , 1960 ; Migdal , 1960 ), even predating 



the 19 67 discovery of pulsars (Hewish et ai, 1968), which were soon identified as rapidly rotating magnetic neutron 
stars (Gold, 1969). Interest in nucleonic pairing has intensified in recent years, owing primarily to experimental 
developments on two different fronts. In the field of astrophysics, a series of X-ray satellites (including Einstein, 
EXOSAT, ROSAT, and ASCA) has brought a flow of data on thermal emission from neutron stars, comprising both 
upper limits and actual flux measurements. The recent launching of the Chandra X-ray observatory provides further 
impetus for more incisive theoretical investigations. On the terrestrial front, the expanding capabilities of radioactive- 
beam and heavy-ion facil ities have stimulated a concerted ex p lorat ion of nuclei far from stability, with a special focus 
on neutron-rich species (Mueller and Sherril, 1993; Riisagei, 1994). Pairing plays a prominent role in modeling the 
structure and behavior of these newly discovered nuclei. 

Since the field is quite vast, we limit our discussion to several recent advances that have taken place. We will focus 
in particular on two overlapping questions: (i) how does many-body pairing evolve from the bare nucleon-nuclcon 
interaction, and (ii) what ar e the experim ental (and perhaps theoretical) manifestations of pairing in finite nuclei? 



195C| ) pointed out that a short-ranged, attractive, nucleon-nucleon interaction 



Over fifty years ago, Mayer (Mayer, 

would yield J = ground states. The realistic bare nucleon-nucleon potential indeed contains short-range attractive 
parts (particularly in the singlet-S* and triplet-P channels) that give rise to pairing in infinite nuclear matter and 
nuclei. In this Review, we will discuss various calculations that demonstrate this effect. We will also demonstrate the 
link between superfiuidity in nuclear matter and its origin from realistic nucleon-nucleon interactions. We then study 
the nature of pair correlations in nuclei, their potential impact on nuclear structure experiments, and the origin of 
pairing in the presence of a random two-body interaction. We conclude with recent experimental evidence that points 
to a relation between pairing and phase transitions (or transformations) in finite nuclear systems. 

Before we present the outline of this work, we feel that some historical remarks about the particularity of the pairing 
problem in nuclear physics may be appropriate. 



A. Theory of pairing in nuclear physics 

In 1911 Kamerlingh Onnes discovered superconductivity in condensed matter systems, and its microscopic expla- 
nation came about th r ough the highly successful pairing theory proposed in 1957 by Bard een, Cooper, and Schrieffer 
(BCS) (ICoopcr et al.[ |l957| ). A series of first applications to nuclear structure followed (|Belyaev|, |l959|; [Bohr et 



1958| ; Migdal r |l959| ). The BCS theory also generalized the seniority coupling scheme in which pair- wise coupling of 



equivalent nucleons to a state of zero angular momentum takes place. The scheme had been developed during years 



previous to the discovery of BCS theory (Mayer, 195C; Racah, 1942; Racah and Talmi, 1953) 



BCS applications in nuclear structure calcu 



ations incorporate two inherent drawbacks. First, the BCS wave 
function is not an eigenstate of the number operator, so that number fluctuation is an issue. Second, there is a critical 
value of the pairing-force strength for which no non-trivial solution exists. Several attem pts were made to overcom e 
these problems: calculating the rando m phase approximation (RPA) in addition to BCS ( Unna and Wenesei , 1965); 
including particle number projection (Kerman et al., 1961) after variation, valid for pairing strengths above the 
critical pairing strength, and a projection before variation that works well for all pairing strength values. A simplified 
prescription for the latter is a technique known as the Lipkin-Nogami method ( Lipkii:i| , 1960; Nogami, 1964). It has 
been quite successful in overcoming some of the shortfalls that occur when BCS is applied to nuclei; see e.g., the recent 
works of Hagino et al. ( plagino and Bertsch , 2000; Hagino et al., 2002) and references therein. Of course, BCS is an 
approximate solution to the many-body problem and assumes a particular form for the many-body wave function. 
Another, more drastic, approximation to the many-body problem assumes that a single Slater determinant suffices to 
describe the nuclear ground state. This mean-field solution to the many-body problem gives rise to Hartree-Fock (HF) 
theory. An effective nucleon-n uclcon potential describes the nuclear interaction, and is typically a parameterization 
of the Skyrme zero-range force (Skyrme, 1956, 1959; Vautherin and Brink, 197C, 1972). Solutions of the HF equations 
describe various nuclear ground-state properties sufficiently (Quentin and Flocard, 1978), but they do not i nclud e 
an explicit pairing interaction. Finite-range interactions, such as the Gogny interaction (Decharge et al., 1975), 



when used in Hartree-Fock calculations has also no pairing by construction. A general way to include pairing into a 
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mean- field description genera ted by e.g., a Skyrme interaction requires solving the Hartree-Fock-Bogoliubov (HFB) 



equations ( Bogolyubov, 195£ ). Recent a pplications to both stable and weakly bound nuclei may be found in, e.g. 



(Dobaczewski et ai, 1996; Duguct et ai, |2002a| Jb|). A renormali zation s c heme for the HFB equation s was recently 
proposed by Bulgac and Yu for a zero range pairing interaction ( Bulgac , 2002 ; Bulgac and Yu , 2002). Rather than 
solving the full HFB equations, one may first calculate the Hartrec - Fock sinele-par ticlc wave functions and use these 
as a basis for solving the BCS equations ( Nayak and Pearson , 1995 ; Tondcui , 1979 ). For stable nuclei with large one- 
or two-neutron separation energies, the HF-fBCS approximation to HFB is valid, but the technique is not able to 
adequately address weakly bound nuclei due to the development of a particle (usually neutron) gas on or near the 
nuclear surface. 

While nuclear mean-field calculations represent a well-founded method to describe nuclear properties, their results do 
not repre sent complete solutio ns to the nuclear many-body problem. Short of a complete solution to the many-nucleon 
problem (Pudliner et al, 1995| ), the interacting shell model is widely regarded as the most broadly capable description 
of low-energy nuclear structure and the one most directly traceable to the fundamental many-body problem. While 
this is a widely accepted statement, applications of the shell model to finite nuclei encounter several difRculties. Chief 
among these is the choice of the interaction. A second problem involves truncations of the Hilbert space, and a third 
problem involves the numerics of solving extremely large eigenvalue problems. 

Skyrme and Gogny forces are parameterized nuclear forces, but they lack a clear link to the bare nucleon-nucleon 
interaction as described by measured scattering phase s hifts. The s ame p hilosophy has been used for shell-model 



interactions, e.g., with the USD Is-Od-shell interaction ( Wildenthal , 1984b ). While quite successful, these types of 
interactions cannot be related directly to the nucleon-nucleon interaction either. The shell model then becomes a true 
model with many parameters. Alternatively, many attempts have been made to derive an effective nucleon-nucleon 
interaction in a given shell-model space from the bare nucleon- nucleon interaction using many-body perturbation 



theory. (For a modern exposition on this difficult problem, see (Hjorth- Jensen et al, 1995) and references therein.) 



While this approach appears to work quite well for many nuclei, there are several indications (Pieper et al, 2001 



Pudliner et al 



1997, 1995 ) that an effective interaction based on a two- body force only fails to reproduce experimental 



data. As shown in e.g., ( Pieper et al , 2001; Pudliner et al , 1997 , 1995), these difficulties are essentially related to the 
absence of a real three-body interaction. It should be noted, however, that the deficiencies of the effective interactions 
are minimal and affect the ground-state energies more than they affect the nuclear spectroscopy. Thus, understanding 
various aspects of physics from realistic two-body interactions, or their slightly modified, yet more phenomenological, 
cousins, is still a reasonable goal. 



B. Outline 



This work starts with an overview of pairing in infinite matter, with an emphasis on superfluidity and superconduc- 
tivity in neutron stars. As an initial theme, we focus on the link between superfluidity in nuclear matter and its origin 
from realistic nucleon-nucleon interactions. This is done in Sec. || where we discuss pairing in neutron star matter 
and symmetric nuclear matter. Thereafter, we focus on various aspects of pairing in finite nuclei, from spectroscopic 
information in Sec. Ill to pairing from random interactions in Sec. |l^ and thermodynamical properties in Sec. 0. 
Concluding remarks are presented in Sec. VI. The paragraphs below serve as an introduction to the exposed physics. 



1. Pairing in neutron stars 



The presence of neutron superfluidity in the crust and the inner part of neutron stars are considered well established 
in the physics of these compact stellar objects. To a first approximation, a neutron star is described as a neutral 
system of nucleons (and possibly heavier baryons) and electrons (and possibly muons) in beta equilibrium a t zero 
tempe ratur e with a central density several tirn e s the s atura t ion dens i ty pp of symmetrical nuclear matter (|Alpar 
|1995| ; IHeiselberg and Hjorth-Jensenj |2000| ; |Lamb| , |l99l| ; |Pcthick| , |l992| ; [Shapiro and Tcukolskyj , |1983| ; [Wiringa 
198§| ). The gross structure of the star (mass, radius, pressure, and density profiles) is determined by solving the 



et a 
et a 

Tolman-Oppenheimer-Volkov general relativistic equation of hydrostatic equilibrium, consistently with the continuity 
equation and the equation of state (which embodies the microscopic physics of the system). The star contains (i) 
an outer crust made up of bare nuclei arranged in a lattice interpenetrated by relativistic electrons, (ii) an inner 
crust where a similar Coulomb lattice of neutron-rich nuclei is embedded in Fermi seas of relativistic electrons and 
neutrons, (iii) a quantum fluid interior of coexisting neutron, proton, and electron fluids, and finally (iv) a core region 
of uncertain constitution and phase (but possibly containing hyperons, a pion or kaon condensate, and/or quark 
matter). Fig. |l] gives a schematic portrait of a possible neutron star structure. 

In the low-density outer part of a neutron star, the neutron superfiuidity is expected mainly in the attractive singlet 
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^Sq channel. Qualitatively, this phenomenon can be understood as follows. At the relatively large average particle 
spacing at the "low" densities involved in this region, i.e., p ~ Po/lO with pq the saturation density of symmetrical 
nuclear matter, the neutrons experience mainly the attractive component of the ^5*0 interaction; however, the pairing 
effect is quenched at higher densities, ~ po and beyond, due to the strong repulsive short-range component of this 
interaction. At higher density, the nuclei in the crust dissolve, and one expects a region consisting of a quantum liquid 
of neutrons and protons in beta equilibrium. By similar reasoning, one thus expects proton pairing to occur 
in the quantum fluid interior, in a density regime where the proton contaminant (necessary for charge balance and 
chemical equilibrium) reaches a partial density pp ~ po/10. In this region, neutron superfluidity is expected to occur 
mainly in the coupled ^P2-^-f2 two-neutron channel. At such densities, one may also expect superfluidity from other 
baryons such as, e.g., hyperons to arise. The possibility for hyperon pairing is an entirely open issue; see, for example, 
(Balberg and Barnea, 1997). Neutron, proton, and eventual hyperon superfluidity in the ^Sq channel, and neutron 
superfluidity in the '^P2 channel, have been shown to occur with gaps of a few MeV or less (Baldo et ai, 1998b); 



however, the density ranges in which gaps occur remain uncertain. In the core of the star any superfluid phase should 
finally disappear, although the possibility of a color superconducting phase may have interesting consequences. At 
large bary on densities f o r whi ch perturbative QCD applies, pairing gaps for like quarks have been estimated to be a 
few MeV (D. and Love, 1984). However, the pairing gaps of unlike qu arks (ud, us, and d s) have been suggested to 



be several tens to hundreds of MeV through no n-perturbative studies ( Alford et al , 1999 ) kindling interest in quark 
superfluidity and superconductivity (Son, 1999| ) and their effects on neutron stars. 

A realistic ab initio prediction of the microscopic physics of nucleonic superfluid components in the interiors of 



neutron stars is crucial to a quantitative understand ing of neutrino cooling mechanisms (Friman and Maxwell, 1979 



Takatsuka and Tamagaki, 1997; Tsuruta, 1979, 1998) that operate immediately after their birth in supernova events, 
as well as the magnetic properties, vortex structure, rotational dynamics, and pulse timing irregularities of these 
superdense stellar objects. In particular, when nucleonic species enter a superfluid state in one or another region of 
the star, suppression factors of the form exp^—Ap/ksT) are introduced into the expression for the emissivity, Ap 
being an appropriate average measure of the energy gap at the Fermi surface. Pairing thus has a major effect on 
the star's thermal evolution through suppressions of neutrino emission processes and specific heats as well; see, for 



example, (Page et al, 200C) 



2. Pairing phenomena in nuclei 



After the excursion to infinite matter, we return to the question concerning how to obtain informatio n o n pairing 
correlations in finite nuclei from abundantly available spectroscopic data. We discuss this point in Sec. III. Even in 
the presence of random interactions, signatures of pairing still remain in finite many-body systems. In Sec. IV we 
present a discussion of pairing derived from random interactions. 

Apart from relatively weak electric forces, the interactions between two protons are very similar to those between two 
neutrons. This yields the idea of charge symmetry of the nuclear forces. Furthermore, the proton-neutron interaction 
is also very similar. This led very early to the idea of isotopic invariance of the nucleon-nucleon interaction. A nucleon 
with quantum number isospin t = 1/2 may be in one of two states, = —1/2 (proton) or — +1/2 (neutron). Of 
course, the symmetry is not exact, but is widely employed when discussing nuclei. It leads to a quantum number T 
called isospin, and its projection ~ {N — Z)/2, where the number of neutrons (protons) in the nucleus is N (Z). 

We can define with this isospin symmetry two distinct states within the two-nucleon system. A T = 1 nucleon- 
nucleon system can have spin-projection = 1, 0, — 1. = 1 corresponds to a neutron- neutron system, = to a 
proton- neutron system, and = — 1 to a proton-proton system. The nucleons in this case have total spin J = in 
order for the full wave function to maintain antisymmetry of the total nucleon-nucleon wave function. For the same 
reason, T = proton- neutron systems can only have = and J = 1. Thus, two different types of elementary 
particle pairs exist in the nucleus, and they depend on both the spin and isospin quantum numbers of the two-particle 
system. 

This brief discussion of the general quantum numbers of a two-nucleon system is a natural starting point for a 
discussion of pairing found in nuclei. All even-even nuclei have a ground-state with total angular momentum quantum 
number and parity, tt, J'^ = 0"*". One can postulate a pairing interaction that couples particles in time-reversed states. 
Using this type of simple pairing interaction, one can also understand the fact that in even-even nuclei the ground 
state is rather well separated from excited states, although in the even-odd neighbor nucleus, several states exist near 
the ground state. 

The behavior of the even-even ground state is usually associated with isovector (T = 1) pairing of the elementary 
two-body system. Simplified models of the nucleon-nucleon interaction, such as the seniority model (Talmi, 1993), 
predict a pair condensate in these systems. An open question concerns evidence for isoscalar (T — 0) pairing in 
nuclei. One unique aspect of nuclei with N=Z is that neutrons and protons occupy the same shell-model orbitals. 
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Consequently, the large spatial overlaps between neutron and proton single-particle wave functions are expected to 
enhance neutron-proton (np) correlations, especially the np pairing. 

At present, it is not clear what the specific experimental fingerprints of the np pairing are, whether the np correla- 
tions are strong enough to form a static condensate, and what their main building blocks are. Most of our knowledge 
about nuclear pairing comes from nuclei with a sizable neutron excess where the isospin T—l neutron- neutron (nn) 
and proton-proton (pp) pairing dominate. Now, for the first time, there is an experimental opportunity to explore 
nuclear systems in the vicinity of the N=Z line which have many valence np pairs; that is, to probe the interplay 
between the like-particle and neutron-proton (T=0,1, Tz=0) pairing channels. One evidence related to T = pairing 
involves the Wigner energy, the extra binding that occurs in N ^ Z nuclei. We will discuss this in more detail in 
Sec. m. 

On e possible way to experimentally access pair correlations in nuclei is by neutron-pair transfer, see e.g., (Yoshida, 
1962| ). Simply stated, if the ground-state of a nucleus is made of BCS pairs of neutrons, then two- neutron transfer 



should be enhanced when compared to one-neut ron transfer . Coll ective enhancement of pair transfer is expected if 
nuclei with ope n shells are brought into contact (Peter et at, 1999). Pairing fluctuations are also expected in rapidly 



rotating nuclei ( Shimizu et al. , 1989 ). I n lighter systems, s uch a s He, two-neutron transfer has been used for studying 



the wave function of the ground state (Oganessian et al, 1999) 

Finally, we discuss phenomenological descriptions of nuclear collective motion where the nuclear ground state and 
its low-lying excitations are represented in terms of bosons. In one such model, the Interacting Boson Model (IBM), 



L = ( S) an d L — 2 (D) bosons are identified with nuclcon pairs having the same quantum numbers (lachello and 



Ariipi], 1988), and the ground state can be viewed as a condensate of such pairs. Shell- model studies of the pair 
structure of the ground state and its variation with the number of valence nucleons can therefore shed light on the 
validity and microscopic foundations of these boson approaches. 



3. Thermodynamic properties of nuclei and level densities 



The theory of pairing in nuclear physics is also strongly related to other fields of physics, such as distinct gaps in 
ultrasmall metallic grains in the solid state. These systems share in common the fact that the energy spectrum of 
a system of part icles confined to a small region is quantized . It is only recently, through a series of experiments by 
Tinkham et al. ( [Black et al , 1996, 1997] ; Ralph et al., 1995), that spectroscopic data on discrete energy levels from 
ultrasmall metallic grains (with sizes of the order of a few nanometers and mean level spacings less than millielectron- 
volts) has been obtained by way of single-electron-tunneling spectroscopy. Measurements in solid state have been 
much more elusive due to the size of the system. The discrete spectrum could not be resolved due to the energy 



scale set by temperature. Of i nterest here is the observation of so-called parity effects. Tinkham et al. (Black et al. 



1996 , 1997 ; Ralph et al. , 1995 ) were able to observe the number parity (odd or even) of a given grain by studying the 



evolution of the discrete spectrum in an applied magnetic field. These effects were also observed in experiments on 
large Al grains. It was noted that an even grain had a distinct spectroscopic gap whereas an odd grain did not. This 
is clear evidence of superconducting pairing correlations in these grains. The spectroscopic gap was driven to zero by 
an applied magnetic field; hence the paramag netic breakdown of pairing correlations could be studied in detail. For 
theoretical interpretations, see, f or example, ( Balian et al. , 1999; von Delft and Ralph , 2001 ; Dukelsky and Sierra , 
1999| ; pVIastellonc and Falci| , |1998| ). 



In the smallest grains with sizes less than 3 nanometers, such distinct spectroscopic gaps could however not be 
observed. This vanishing gap revived an old issue: what is the lower size of a system for the existence of supercon- 
ductivity is such small grains? 

A nucleus is also a small quantal system, with discrete spectra and strong pairing correlations. However, whereas the 
statistical physics of the above experiments on ultrasmall grains can be well described through a canonical ensemble, 
i.e., a system in contact with a heat bath, the nucleus in the laboratory is an isolated system wi th no heat exchang e 
with the environment. The appropriate ensemble for its description is the microcanonical one (Balian et al, 1999). 
This poses significant interpretation problems. For example, is it possible to define a phase transition in an isolated 
quantal system such as a nucleus? 

In Sec. we attempt to link our discussion to such topics via recent experimental evidence of pairing from studies 
of level densities in rare-earth nuclei. The nuclear level density, the density of eigenstates of a nucleus at a given 
excitation energy, is the important quantity that may be used to describe thermodynamic properties of nuclei, such 
as the nuclear entropy, specific h eat, and temp erature. Bethe first described the level density using a non-interacting 
fermi gas model for the nucleo ns (Bcthc , 1936). Mod i fications to th i s pict ure, such as the back-shifted fermi gas which 
includes pair and shell effects ( Gilbert and Cameron| , 1965; Newton, 1956| ) not present in Bethe's original formulation, 
are in wide use. These modifications incorporate long-range pair correlations that play an important role in the low 
excitation region. Experimentalists recently developed methods (Henden et al, 1995; Tveter et al, 1996) to extract 
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level densities at low spin from measured 7-spectra. 

There is evidence for the existence of paired nucleons (Cooper pairs) at low temperature^. In high-spin nuclear 
physics, the backbending phenomenon is a beautiful manifestation of the breaking of pairs. The mechanism induced 



by Coriolis forces tends to align single partic l e angular momenta a l ong th e nuclear rotational axis ( Faessler et al 



1976 



Johnson et al, 1971; Ricdingcr et al, 1980; Stephens and Simon, 1972). Theoretical models also predict a reduction 



in th e pair correlations at higher temperatures (|D0ssing et al , 1995 ; Mottelson and Valantin , 1960 ; Muhlhans et 



1983). There is also an interesting connection between quasiparticle spectra in metallic grains and high-spin spectra 
in nuclei. In nuclei it is the Coriolis force that acts on pairs of nucleons and plays thus a role similar to the magnetic 
field acting on Cooper pairs of electrons. 



The breaking of pairs is difficult to observe as a function of intrinsic excitation energy. Recent theoretical (D0ssing 
et ) and experimental ( [Melby et al , 1999; Tvetcr et al, 1996| ) works indicate that the process of breaking 



pairs takes place over several MeV of excitation energy. Thus, the phenomenon of pair breaking in a finite-fermi 
system behaves somewhat differently than what would be expected in nuclear matter. The corresponding critical 
temperature in finite systems is measured to be Tc ~ 0.5 MeV/fcs ( ^chiller et al , 2001), where ks is Boltzmann's 
constant. Recent work extracted the entropy of the I6i,i62-|-)y ^^^^ iYi.iy^Yb isotopes and deduced the number of 
excited quasiparticles as a function of excitation energy. We describe this result in more detail in Sec. 0. 



II. PAIRING IN INFINITE MATTER AND THE NUCLEON-NUCLEON INTERACTION 

Pairing correlations and the phenomenon of superconductivity and superfluidity are intimately related to the un- 
derlying interaction whether it is, for example, the nucleon-nucleon (NN) interaction or the interaction between '^He 
atoms. In this section we discuss, through simple examples, some of the connections between pairing correlations as 
they arise in nuclear systems and the bare NN interaction itself, that is, the interaction of a pair of nucleons in free 
space. The latter is most conveniently expressed in terms of partial waves (and their pertaining quantum numbers 
such as orbital angular momentum and total spin) and phase shifts resulting from nucleon-nucleon scattering exper- 
iments. Actually, without specializing to some given fermionic systems and interactions, it is possible to relate the 
pairing gap and the BCS theory of pairing to the experimental phase shifts. This means, in turn, that we can, through 
an inspection of experimental scattering data, understand which partial waves may yield a positive pairing gap and 



eventually lead to, e.g., a superfluid phase transition in an infinite fermionic system. We show this in subsec. [I.C 
(although we limit the attention to nuclear interactions), after we have singled out those partial waves and interaction 
properties which are expected to be crucial for pairing correlations in both nuclei and neutron stars. These selected 
features of the NN interaction are discussed in the next subsection. A brief overview of superfluidity in neutron stars 



and pairing in symmetric nuclear matter is presented in subsec. II. D, with an emphasis on those partial waves of 



the NN interaction which are expected to produce a finite pairing gap. Features of neutron-proton pairing in infinite 



matter are reviewed in subsec. II. D. 2. Concluding remarks, open problems, and perspectives are presented in the last 
subsection. 



A. Selected features of the nucleon-nucleon interaction 

The interaction between nucleons is characterized by the existence of a strongly repulsive core at short distances, 
with a characteristic radius ^ 0.5 — 1 fm. The interaction obeys several fundamental symmetries such as translational, 
rotational, spatial-reflection, time-reversal invariance and exchange symmetry. It also has a strong dependence on 
quantum numbers such as total spin S and isospin T, and, through the nuclear tensor force which arises from, e.g., 
one-pion exchange, it also depends on the angles between the nuclcon spins and separation vector. The tensor force 
thus mixes different angular momenta L of the two-body system, that is, it couples two-body states with total angular 
momentum J = L — 1 and J = L + 1. For example, for a proton-neutron two-body state, the tensor force couples the 
states ■^S'l and ^Di, where we have used the standard spectroscopic notation '^^'^^Lj. 

Although there is no unique prescription for how to construct an NN interaction, a description of the interaction 



in terms of various meson exchanges is presently the most quantitative representation, see for example (Machleidt, 
1989| , |2001| ; [Machleidt et al\ , [1996[ ; ^toks et al\ , [1994[ ; [Wiringa et a/] , [1995| ), in the energy regime of nuclear structure 



The concept of temperature in a microcanonical system such as the nucleus is highly non-trivial. Temperature itself is defined by a 
measurement process, involving thereby the exchange of energy, a fact which is in conflict with the definition of the microcanonical 
ensemble. It is only in the thermodynamic limit that e.g., the caloric curves in the canonical and microcanonical ensembles agree. The 
word temperature in nuclear physics should therefore be used with great care. 
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physics. We will assume that meson-exchange is an appropriate picture at low and intermediate energies. Further, in 
our discussion of pairing, it suffices at the present stage to limit our attention to the time-honored configuration-space 
version of the nucleon-nucleon (NN) interaction, including only central, spin-spin, tensor and spin-orbit terms. In 
our notation below, the mass of the nucleon Mat is given by the average of the proton and neutron masses. The 
interaction reads (omitting isospin) 



F(r) ^\C^ + Ch + C^ai • era + Ct 1 



{mar) 



1 



1 



SL 



(niary 



L S 



niaT 



where ma is the mass of the relevant meson and 5*12 is the tensor term 



(1) 



(2) 



where cr is the standard operator for spin 1/2 particles. Within meson-exchange models, we may have, e.g., the 
exchange of tt, ry, p, cr, and 5 mesons. As an example, the coefficients for the exchange of a tt meson are — Ct — 

2 ^ 
SnN-k ^ 



4ir 12Mi 



and 



CsL — with the experimental value for g%j\j^ sa 13 — 14; see, for example, ( Machleidt , 



2001 ) for a recent discussion. 

The pairing gap is determined by the attractive part of the NN interaction. In the ^5*0 channel the potential is 
attractive for momenta k < 1.74 fm^^ (or for interparticle distances r > 0.6 fm), as can be seen from Fig. ^. In the 
weak coupling regime, where the interaction is weak and attractive, a gas of ferniions may undergo a superconducting 
(or superfluid) instability at low temperatures, and a gas of Cooper pairs is formed. This gas of Cooper pairs will be 
surrounded by unpaired fermions and the typical coherence length is large compared with the interparticle spacing, 
and the bound pairs overlap. With weak coupling we mean a regime where the coherence length is larger than the 
interparticle spacing. In the strong-coupling limit, the formed bound pairs have only a small overlap, the coherence 
length is small, and the bound pairs can be treated as a gas of point bosons. One expects then the system to undergo 
a Bose- Einstein condensation into a single quantum state with total momentum fc = (Nozieres and Schmitt-Rink , 
1985). For the channel in nuclear physics, we may actually expect to have two weak-coupling limits, namely 



when the potential is weak and attractive for large interparticle spacings and when the potential becomes repulsive 
at r « 0.6 fm. In these regimes, the potential has values of typically some few MeV. One may also loosely speak 
of a strong-coupling limit where the NN potential is large and attractive. This takes place where the NN potential 
reaches its maximum, with an absolute value of typically ~ 100 MeV, at roughly ~ 1 fm, see again Fig. |^. We note 
that fermion pairs in the ^5*0 wave in neutron and nuclear matter will not undergo the above-mentioned Bose-Einstein 
condensation, since, even though the NN potential is large and attractive for certain Fermi momen ta, the coherence 
length will always be larger than the interparticle spacing, as demonstrated by De Blasio et al. (De Blasio et al 



1997). The inclusion of in-medium effects, such as screening terms, are expected to further reduce the pairing gap and 



thereby enhance further the coherence length. This does not imply that such a transition is not possible in nuclear 
matter. A recent analysis by Lombardo et al. (Lombardo et al, 2001a; Lombardo and Schuckj , |2001| ) of triplet ^Si 
pairing in low-density symmetric and asymmetric nuclear matter indicates that such a transition is indeed possible. 

Hitherto we have limited our attention to one single partial wave, the ^Sq channel. Our discussion about the 
relation among the NN interaction, its pertinent phase shifts, and the pairing gap, can be extended to higher partial 
waves as well. An inspection of the experimental phase shifts for waves with J < 2 and total isospin T — 1, 
see Fig. |^, reveals that there are several partial waves which exhibit attractive (positive phase shifts) contributions 
to the NN interaction. Such attractive terms are in turn expected to yield a possible positive pairing gap. This 
means that the energy dependence of the nucleon-nucleon (NN) phase shifts in different partial waves offers some 
guidance in judging what nucleonic pair-condensed states are possible or likely in different regions of a neutron star. 
A rough correspondence between baryon density and A'^A'^ bombardment energies can be established through the 
Fermi momenta assigned to the nucleonic components of neutron-star matter. The lab energy relates to the Fermi 
energy through E^i, = itp — 4?i^fc|,/2Mjv. This is demonstrated in Fig. ^ for various NN interaction models that fit 
s cattering data up to E 'lab ~ 350 MeV. For c o mpar ison, we include results for older potential models such as the Paris 
( Lacombe et al , 198C), V14 ( Wiringa et al. , 1984 ) and Bonn B ( Machleidt , 1989) interactions. Note, as well, that 
beyond the point where these potential models have been fit, there is a considerable variation. This has important 
consequences for reliable predictions of the '^P2 pairing gap. 

In pure neutron matter, only T = 1 partial waves are allowed. Moreover, one need only consider partial waves with 
L < 4 in the range of baryon density - optimistically, p < (3 — A)pQ- where a nucleonic model of neutron-star material 
is tenable, where po = 0.16 fm~'^ is the saturation density of nuclear matter. We have already seen that the ^Sq phase 
shift is positive at low energy (indicating an attractive in-medium force) but turns negative (repulsive) at around 250 
MeV lab energy. Thus, unless the in-medium pairing force is dramatically different from its vacuum counterpart, the 
situation already suggested above should prevail: S-wave pairs should form at low densities but should be inhibited 
from forming when the density approaches that of ordinary nuclear matter. 
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The next lowest T = 1 partial waves are the three triplet P waves ^Pj, with J = 0,1,2. For the the ^Pq state, 
the phase shift is positive at low energy, turning negative at a lab energy of 200 MeV. The attraction is, however, 
not sufficient to produce a finite pairing gap in neutron star matter. The ^Pi phase shift is negative at all energies, 
indicating a repulsive interaction. The phase shift is positive for energies up to 1 GeV and is the most attractive 
T = 1 phase shift at energies above about 160 MeV. Whereas the ^5*0 partial wave is dominated by the central force 
contribution of the NN interaction, see Eq. (|l|), the main contribution to the attraction seen in the partial wave 
stems from the two-body spin-orbit force for intermediate ranges in Eq. (|l]), i.e., the term proportional with L ■ S. This 
is demonstrated in Fig. ^ where we plot the coordinate space version of the Argonne Vi^ interaction ( Wiringa et al. 



1995) with and without the spin-orbit contribution. Moreover, there is an additional enhancement due to the ^P2- P2 



tensor force. A substantial pairing effect in the P2- P2 channel may hence be expected at densities somewhat in 
excess of p^, again assuming that the relevant in- vacuum interaction is not greatly altered within the medium. 

The remaining T = 1 partial waves with L < 4 are both singlets: and ^G^. However, the phase shifts of these 
partial waves, albeit being positive over the energy domain of interest, do not provide any substantial contribution to 
the pairing gap. Thus, only the and '^P2 partial waves yield enough attraction to produce a finite pairing gap in 
pure neutron matter. Singlet and triplet pairing are hence synonomous with and ^P2-'^P2 pairing, respectively. 



B. Pairing gap equations 

The gap equation for pairing in non-isotropic partial waves is, in general, more complex than in the simplest singlet 
iS-wave case, in particular in neutron and nuclear matter, where the tensor interaction can couple two different partial 
waves (Baldo et al, 1995; Takatsuka and Tamagaki, 1993; Tamagaki, 1970| ). This is indeed the situation for the ■^P2- 
^P2 neutron channel or the '^Si-'^Di channel for symmetric nuclear matter. For the sake of simplicity, we disregard 
for th e moment spin degrees of freedom and the tensor interaction. Starting with the Gorkov equations (^chrieffei , 



1964), which involve the propagator G(fc,a;), the anomalous propagator F{k,uj), and the gap function A(fc), we have 



c-^fe) -A(fe)\/G\ /1\ 
-At(fc) c^ + e(fc) j l^Ft j ^'"''^^^ 1^0 j ' 

where e(fc) = e{k) — fi, fi being the chemical potential and e{k) the single-particle spectrum. The quasi-particle energy 
E{k) is the solution of the corresponding secular equation and is given by 

i?(fe)2 = 6(fc)2-HiA(fc)r (4) 

The anisotropic gap function A(A;) is to be determined from the gap equation 

Aik) = -Y^ik\V\k')^ (5) 
fc' ^ ^ 

The angle-dependent energy denominator in this equation prevents a straightforward separation into the different 
partial wave components by expanding the potential. 



(k\V\k')^A7Ty^{2L + l)PL{k-k')VL{k,k'), (6) 



L 



and the gap function, 



^C^) - E \l:^YLM{k)ALM{k), (7) 



. 21 + 1 

LM 



with L and M being the total orbital momentum and its projection, respectively. The functions Yl,m are the spherical 
harmonics. However, after performing an angle average approximation for the gap in the quasi-particle energy, 

|A(fc)p D{kf ^^jdk |A(fc)p = _L_|AiM(A;)P , (8) 

the kernels of the coupled integral equations become isotropic, and one can see that the different m-components 
become uncoupled and all equal. One obtains the following equations for the partial wave components of the gap 
function: 

Ai(fc) = -- r k'dk' ^ VLik,k') 
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Note that there is no dependence on the quantum number M in these equations; however, they still couple the 
components of the gap function with different orbital momenta L (^5*0, ^Pqj ^Pi^ '^P2, ^^2, ^^2, etc. in neutron 
matter) via the energy denominator. Fortunately, in practice the different components Vl of the potential act mainly 
in non-overlapping intervals in density, and therefore also this coupling can usually be disregarded. 

The addition of spi n degrees of freedom and of the tensor force does not change the picture qualitatively and is 



explained in detail in ( Baldo et ai , 1995 ; Takatsuka and Tamagaki 



1993|) . The only modification is the introduction 
of an additional 2x2 matrix structure due to the tensor coupling of the ■^P2 and channels. Such coupled channel 
equations can be written as 



(fc) 

Eikf 
D{kf 



1 



dk'k 



1 



Vll 

[e{k) - e{kp)f + D{k)\ 



-Vll' 
Vl'u 



{k,k') 



Al 



ik'), 



(10) 

(11) 
(12) 



Here e(fc) — k^ /2m + U{k) arc the single-particle energies of a neutron with momentum k, and kp is the Fermi 
momentum. The orbital momenta L and L' could, e.g., represent the '^P2 and ^^^2 channel, respectively. Restricting 
the attention to only one partial wave, it is easy to get the equation for an uncoupled channel like the ^Sq wave, i.e., 
we obtain 



1 r°° 

A(fc)L = / dk'k'^VLL{k,k' 

I" Jo 



,^A(fc') 



E{k' 



(13) 



where Vphik^k') is now the bare momentum-space NN interaction in the ^So channel, and E(k) is the quasiparticle 
energy given by E{k) — (e(fc) — e{kp))^ + A(fc)|. 
The quantities 



Vll' ik,k')= / drr^jL, ik'r)VLL' {r)jL (kr) 



(14) 



are the matrix elements of the bare interaction in the different coupled channels, e.g., (T = 1; S ^ 1; J — 2; L, L' — 
1,3). It has been shown that the angle average approximat i on is an excellent approximation to the true solution that 
involves a gap function with ten components (Kodel et ai, 1996; Takatsuka and Tamagaki] , 1993), as long as one is 
only interested in the average value of the gap at the Fermi surface, A^ = D(kp), and not the angular dependence 
of the gap functions Ai(A;) and Ai/(fc). 

Recently Kho del, Khodel, and Cl ark (Kodel et ai, 1998| , |2001| ) proposed a separation method for the triplet pairing 
gap, based on (Kodel et al, 1996[ ), which allows a generalized solution of the BCS equati o n tha t is numeric ally 



reliable, without employing an angle-average approach. We refer the reader to ( Kodel et al , 199(j , 1998 , 2001 ) for 
more details. In this approach, the pairing matrix elements are written as a separable part plus a remainder that 
vanishes when either momentum variable is on the Fermi surface. This decomposition effects a separation of the 
problem of determining the dependence of the gap components in a spin-angle representation on the magnitude of 
the momentum (described by a set of functions independent of magnetic quantum number) from the problem of 
determining the dependence of the ga p on angle o r magnetic p rojection. The former problem is solved through a set 



of nonsingular, quasilinear equations ( Kodel et al , 1998 , ^001 ). There is, in general, a good agreement between their 
approach and the angle average scheme. However, the general scheme of Khodel, Khodel, and Clark offers a much 
more stable algorithm for solving the pairing gap equations for any channel and starting with the bare interaction 
itself. In nuclear physics the interaction typically has a strongly repulsive core, a fact that can complicate significantly 
the iterative solution of the BCS equations. 

An important ingredient in the calculation of the pairing gap is the single-particle potential U{k). The gap equation 
is extremely sensitive to both many-body renormalizations of the interaction and the similar corrections to the single- 
particle energies. Many-body renormalizations of the interaction will be discussed in Sec. [I.E[ In our discussion 
below, we will present results for various many-body approaches to U{k), from U(k) = to results with different 
Brueckner-Hartree-Fock (BHF) calculat ions, with both a dis continuous choice, a model-space BHF approach and 
within the "continuous-choice" scheme (Icukennc et al, 1976| ). 

The single-particle energies appearing in the quasiparticle energies (^ and (|l^) are typically obtained through a 
self-consistent BHF calculation, using a G-matrix defined through the Bethc-Brueckner-Goldstone equation as 



G = V + V 



Lo - Ho 



G, 



(15) 
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where V is the nucleon-nucleon potential, Q is the Pauh operator which prevents scattering into intermediate states 
prohibited by the Pauli principle, Hq is the unperturbed Hamiltonian acting on the intermediate states, and to is the 
s tarting energy, the unpert urbed energy of the interacting particles. Methods to solve this equation are reviewed in 



(Hjorth- Jensen et al, 1995). The single-particle energy for state ki {i encompasses all relevant quantum numbers like 



momentum, isospin projection, spin, etc.) in nuclear matter is assumed to have the simple quadratic form 



2M 



(16) 



N 



where is the effective mass. The terms M'^ and 5, the latter being an effective single-particle potential related to 
the G-matrix, are obtained through the self-consistent BHF pro cedure. The model-space B HF (MBHF) method for 
the single-particle spectrum has also been used, see, for example, ( Hjorth- Jensen et a/. , 1995 ), with a cutoff momentum 



3.0 fm > kp. In this approach the single-particle spectrum is defined by 



with the single-particle potential Ui given by 



k,,<ki 



{kikh\G{uj 
0,k, 



2Mn 



^ ffe, + '^kJlhkh) AS 



(17) 



(18) 



where the subscript AS denotes antisymmetrized matrix elements. This prescription reduces the discontinuity in the 
single-particle spectrum as compared with the standard BHF choice /cm = kp. The self-consistency scheme consists 
of choosing adequate initial values of the effective mass and 5. The obtained G-matrix is then used to calculate the 
single-particle potential m, from which we obtain new values for m* and 5. This procedure continues until these 
parameters vary little. 



Recently, Lombardo et al. (Lombardo et al., 2001b; Lombardo and Schulze, 2001) have reanalyzed the importance 



of the various approaches to the single-particle energies. Especially, they demonstrate that the energy dependence of 
the self-energy can deeply affect the magnitude of the energy gap in a strongly correlated Fermi system; see also the 
recent works of Bozek in (Bozek, 1999, 2000, 2002). We will discuss these effects in Subsec. H.E. 



C. Simple relations between the interaction and the pairing gap for identical particles 

1. The low density limit 



A general two-body Hamiltonian can be written in the form H — Hi + H2 where 

a 



(19) 
(20) 



where and a are fermion creation and annihilation operators, and V are the uncoupled matrix elements of the 
two-body interaction. The sums run over all possible single-particle quantum numbers. 

We limit the discussion in this section to a Fermi gas model with two-fold degeneracy and a pairing-type interaction 
as an example; i.e., the degeneracy of the single-particle levels is set to 2s -I- 1 = 2, with s = 1/2 being the spin of the 
particle. We specialize to a singlet two-body interaction with quantum numbers I = and 5 = 0, that is a ""^ So state, 
with I the relative orbital momentum and S the total spin. For this partial wave, the NN interaction is dominated 
by the central component in Eq. (^, which, within a meson-exchange picture, can be portrayed through 27r (leading 
to an effective a meson) and higher tt correlations in order to yield enough attraction at intermediate distances. 

At low densities, the interaction can be characterized by its scattering length only in order to get expansions 
for the energy density or the excitation spectrum. For the nucleon-nucleon interaction, the scattering length is 
oq = —18.8 ± 0.3 fm for neutron- neutron scattering in the ^Sq channel. If we first assume discrete single-particle 
ener gies , the scattering length approximation leads to the following approximation of the two-body Hamiltonian of 
Eq. @ 



H 



E 



1 



:G 



ij>0 



(21) 
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The indices i and j run over the number of levels L, and the label i stands for a time-reversed state. The parameter G 
is now the strength of the pairing force, while is the single-particle energy of level i. Introducing the pair-creation 
operator = a|m'^i-m' '^^^ rewrite the Hamiltonian in Eq. as 



(22) 



ij>0 



where Ni = aja^ is the number operator, and = id so that the single-particle orbitals are equally spaced at intervals 
d. The latter commutes with the Hamiltonian H . In this model, quantum numbers like seniority S are good quantum 
numbers, and the eigenvalue problem can be rewritten in terms of blocks with good s eniori ty. Loosely speaking, the 
seniority quantum number S is equal to the number of unpaired particles; see (Talmi, 1993| ) for fur ther details . As it 
stands E q. (|2l|) , l ends its elf for shell-model studies. Furthermore, in a series of papers, Richardson ( Richardson , 1963 



19654 ^, 1966a ^, 19674 |b|) obtained the exact solution of the pairing Hamiltonian, with semi- analytic (since there is 



still the need for a numerical solution) expressions for the eigenvalues and eigenvectors. The exact solutions have had 
important consequences for several fields, from Bose condensates to nuclear superconductivity. 

We will come back to this model in our discussion of level densities and thermodynamical features of the pairing 
Hamiltonian in finite systems in Sec. 

Here we are interested in features of infinite matter with identical particles, and using 
rewrite Eq. (^iT) as 



V^/(2^)3/^d3fc, we 



d^k 

W) 



^ekaal^aka + GV 



d^k r d^k 



(2n) 



(27r) 



_a-k'-ak'+- 



(23) 



The first term represents the kinetic energy, with eka — k'^/2m. The label a = ±1/2 stands for the spin, while V is 
the volume. The second term is the expectation value of the two-body interaction with a constant interaction strength 
G. The energy gap in infinite matter is obtained by solving the BCS equation for the gap function A(A;). For our 
simple model we see that Eq. (|l3|) reduces to 



GV 



dk'k'- 



1 



E{k'y 



(24) 



with E{k) the quasiparticle energy given by E{k) — y^{e{k) — e(fci?))2 + A{k)'^, where e(fc) is the single-particle energy 
of a neutron with momentum fc, and kp is the Fermi momentum. Medium effects should be included in e(fc), but we 
will use free single-particle e nergies e(fc) — fc^/2Mjv. 



Papenbrock and Bertsch (Papenbrock and Bertsch, 1999| ) obtained an analytic expression for the pairing gap in 
the low-density limit by combining Eq. (|24|) with the equation for the scattering length oq and its relation to the 
interaction 



MnGV 



+ 1 



GV 



d^k- 



^{e{k)-e{kF)f 



(25) 



which is divergent. However, the authors of ( Papenbrock and Bertsch , 1999| ) showed that by subtracting Eq. (|2j) and 
Eq. (Eq), one obtains 



MnG 
47ran 



G 



2(27r)3 



d^k 



E{k) ^(e(fc)-e(fc^))2 



(26) 



which is no longer divergent. Moreover, we can divide out the interaction strength and obtain 



M 



N 



2(27r) 



d^k 



E{k) ^{e(k)-e{kF)Y 



(27) 



Using dimensional regularization techniques, Papenbrock and Bertsch (Papenbrock and Bertsch, 1999) obtained the 
analytic expression 



^ = (l + X^)V4p,/,(_ 



kpCiQ 



IhJl+X' 



(28) 
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where x = ^{kp) / e{kp) and P1/2 denotes a Legendre function. With a given fermi momentum, we can thus obtain 



the pairing gap. F or smaU values of kpuo, one obtains the weU- known result ( Gorkov and Mclik-Barkhudarov , 1961 
Kodel et al\ , |1996| ) 



A{kp) = -j-Aexp 



2kp\ao\ 



(29) 



This comes about by the behavior of Pi^2{z)7 which has a logarithmic singularity ai z — ~1 (see (Erdelyi, 1953)). For 
large values of kpag, the gap is proportional to e{kp), approaching A « l.l6e{kp). The large val" 
length (flQ = —18.8 ± 0.3fm) clearly limits the domain of validity of the Hamiltonian in Eq. (pT|). 



ue of the scattering 
However, Eq. (|2|) 

provides us with a useful low-density result to compare with results arising from numerical solutions of the pairing 
gap equation. The usefulness of Eq. (|2^) cannot be underestimated: one experimental parameter, the scattering 
length, allows us to make quantitative statements about pairing at low densities. Polarization effects arising from 
renormalizations of the in-medium effective interaction can however change this behavior, as demonstrated recently 
in (Heiselberg et ai, 2000; 3chulze et al, 2001) (see the discussion in Subsec. [I.E). 



2. Relation to phase shifts 



With the results from Eq. ( p9| ) in mind, we ask the question whether we can obtain information about the pairing 
gap at higher densities, without resorting to a detailed model for the NN interaction. 

Here we show that this is indeed the case. Through the experimental phase shifts, we show that one can determine 
fairly accurately the ^Sq pairing gap in pure neutron matter without needing an explicit model for the NN interaction . 
It ou ght to be mentioned that this was demonstrated long ago by e.g., Emery and Scssler, see ( Emery and Sesslei , 
196C). Their approach is however slightly different from ours 



As we saw in the previous subsection, a characteristic feature of ^5*0 NN scattering is the large, negative scattering 
length, indicating the presence of a nearly bound state at zero scattering energy. Near a bound state, where the 
NN T-matrix has a pole, it can be written in separable form, and this implies that the NN interaction itself to a 
good approximation is rank-one separable near this pole (Kodel et ai, 1996; Kwong and Kohler, 1997). Thus, at low 
energies, we approximate 



V{k,k') = Xv{k)v{k'), 

where A is a constant. Then it is easily seen from Eq. (O) that the gap function can be rewritten as 



1 = 



dk'k 



E{k') ' 



(30) 



(31) 



Numerically, the integral on the right-hand side of this equation depends very weakly on the momentum structure 
of A(/c), so in our calculations we could take A(fc) « A^ in E{k). Then Eq. (31) shows that the energy gap Ap is 
determined by the diagonal elements Aw^ {k) of the NN interaction. The crucial point is that in scattering theory it can 
be shown that the inverse scattering problem, that is, the determination of a two-particle potential from t he knowledge 
of the phase shifts at all energies, is exactly, and uniquely, solvable for rank-one separable potentials (Chadan and 
Sabatier, 1992). Following the notation of (Brown and Jackson, 1976| ), we have 



Xv^{k) 



+ K% sinJ(fc) 



/fc2 



(32) 



for an attractive potential with a bound state at energy E 



shift as a function of momentum fc, while a{k) is given by a principal value integral: 



In our case w 0. Here 6{k) is the -^Sq phase 



1 f+°° 
a{k) = -P / dk 
7^ J-00 



, S(k') 



(33) 



where the phase shifts are extended to negative momenta through (5(— fc) — —S{k) ( Kwong and Kohlei , 1997). 

From this discussion we see that Xv'^{k), and therefore also the energy gap A/^, is completely determined by the ^So 
phase shifts. However, there are two obvious limitations on the practical validity of this statement. First of all, the 
separable approximation can only be expected to be good at low energies, near the pole in the T-matrix. Secondly, 
we see from Eq. (33) that knowledge of the phase shifts 6{k) at all energies is required. This is, of course, impossible, 
and most phase shift analyses stop at a laboratory energy E'lab = 350 MeV. The ^Sq phase shift changes sign from 
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positive to negative at i5iab ~ 248.5 MeV; however, at low values of kp, knowledge of v{k) up to this value of k may 
actually be enough to determine the value of Ap, as the integrand in Eq. (|3|) is strongly peaked around kp. 

The in put in our calcula tion is the ^5*0 phase shifts taken from the recent Nijmegen nucleon-nucleon pha se shift 



ana lysis (Stoks et al , 1993 ). We then evaluated Xv'^{k) from Eqs. (32) and (|33|), using methods described in ( Brown 



and Jackson, 1976) to evaluate the principal value integral in Eq. (33). Finally, we evaluated the energy gap for 
various values of kp by solving Eq. (^ij), which is an algebraic equation due to the approximation A{k) fa Ap in the 
energy denominator. 

The resulting energy gap obtained from the experimental phase shifts only is plotted in Fig. ^. In the same figure 
we also report the results (dot-dashed line) obtained using the effective range approximation to the phase shifts: 



k cot S{k) = h -rofc^ 



ao 



(34) 



where ao = —18.8 ± 0.3 fm and tq = 2.75 ± 0.11 fm are the singlet neutron-neutron scatterin g length and effective 
range , respectively. In this case an analytic expression can be obtained for Aw^(fc), as shown in ( phadan and Sabatiei , 
199^ ): 



Xv^k) 



Ik^+P, 



k-^-pl 



(35) 



with — —2/arQ, (3i « —0.0498 fm""'^, and 02 ~ 0.777 fm^^. The phase shifts using this approximation are positive 
at all energies, and this is reflected in Eq. ( |35| ) where Xv^{k) is attractive for all k. From Fig. || we see that below 
kp = 0.5 fm~^ the energy gap can, with reasonable accuracy, be calculated with the interaction obtained directly 
from the effective range approximation. One can therefore say that at densities below kp = 0.5 fm~^, and at the 
crudest level of sophistication in many-body theory, the superfluid properties of neutron matter are determined by 
just two parameters, namely the free-space scattering length and effective range. At such densities, more complicated 
many-body terms are also less important. Also interesting is the fact that the phase shifts predict the position of the 
first zero of A{k) in momentum space, since we see from Eq. ( p5| ) that A(fc) = Apv{k) — first for 6{k) = 0, which 
occurs at Ei^h ~ 248.5 MeV (pp sca t tering ) cor responding to fc » 1.74 fm""'^. This is in good agreement with the 
results of Khodel et al. (Kodel et al, 1996| ). In (Kodcl et al, 1996), it is also shown that this first zero of the gap 
function determines the Fermi momentum at which Ap ~ 0. Our results therefore indicate that this Fermi momentum 
is in fact given by the energy at which the ^Sq phase shifts become negative. 

In Fig. |6| we show also results obtained with recent NN interaction models parametrized to reproduce the Nijmegen 
phase shift da ta. We have here employed the CD-Bonn potential (Machleidt et ai, 1996| ), the Nijmegen I and Nimegen 
II potentials (Stoks et al, 1994). The results are virtually identical, with the maximum value of the gap varying from 
2.98 MeV for the Nijmegen I potential to 3.05 MeV for the Nijmegen II potential. As the reader can see, the agreement 
between the direct calculation from the phase shifts and the CD-Bonn and Nijmegen calculation of Ai? is satisfying, 
even at densities as high as kp = l.A fm~^. The energy gap is to a remarkable extent determined by the available ^Sq 
phase shifts. Thus, the quantitative features of ^Sq pairing in neutron matter can be obtained directly from the ^Sq 
phase shifts. This happens because the NN interaction is very nearly rank-one separable in this channel due to the 
presence of a bound state at zero energy, even for densities as high as as kp = 1.4 fm~^ ^. This explains why all bare 
NN interactions give nearly identical results for the ^5*0 energy gap in lowest-order BCS calculations. Combined with 
Eq. (p^, we have a first approximation to the pairing gap with experimental inputs only, phase shifts, and scattering 
length. 

However, it should be mentioned that this agreement is not likely to survive in a more refined calculation, for 
instance, if one includes the density and spin-density fluctuations in the effective pairing interaction or renormalized 
single-particle energies. Other partial waves will then be involved, and the simple arguments employed here will, of 
course, no longer apply. 



^ This is essentially due to the fact that the integrand in the gap equation is strongly peaked around the diagonal matrix elements. 
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D. Superfluidity in neutron star matter and nuclear matter 



1. Superfluidity in neutron star matter 

As we have seen, the presence of two different superfluid regimes is suggested by the known trend of the nucleon- 
nucleon (NN) phase shifts in each scattering channel. In both the ^5*0 and '^P2-^F2 channels the phase shifts indicate 
that the NN interaction is attractive. In particular for the ^Sq channel, the occurrence of the well-known virtual state 
in the neutron-neutron channel strongly suggests the possibility of a pairing condensate at low density, while for the 
^P2-^F2 channel the interaction becomes strongly attractive only at higher energy, which therefore suggests a possible 
pairing condensate in this channel at higher densities. In recent years, the BCS gap equation has been solved with 
realistic interactions, and the results confirm these expectations. 

The ^Sa neutron superfluid is relevant for phenomena that can occur in the inner crust of neutron s tars, l i ke th e 
formation of glitches, which may be related to vortex pinning of the superfluid phase in the solid crust ( ^auls , 1989). 
The results of different groups are in close agreement on the ^S{) pairing gap values and on its density dependence 



whi ch shows a peak value of about 3 MeV at a Fermi momentu m close to fc^ ~ 0.8 fm ( Baldo et al , 1990 ; Elgar0y 



and Hjorth- Jensen, 199^ ; Kodcl et al, 1996; Schulze et al., 1996| ). All these calculations adopt the bare NN interaction 



as the pairing force, and it has been pointed out that the screening by the medium of the interaction could strongly 
reduce the pairing strength in this channel (Ainsworth et al., 1989, 1993; Chen et al., 1986; Schulze et al., 1996| ). The 
issue of the many-body calculation of the pairing effective interaction is a complex one and still far from a satisfactory 



solution (see also the discussion in Sec. 11. E) 



The precise knowledge of the P2- F2 pairing gap is of paramount relevance for, e.g., the cooling of neutron stars, 
and different values correspond to drastically different scenarios for the cooling process. Generally, the gap suppresses 
the cooling by a factor ~ exp(— A/T) (where A is the energy gap), which is severe for temperatures well below the 
gap energy. Unfortunately, only few and partly contra dictory calculations of the '^P2-^F2 pairing g a p exi s t in the 



literat ure, even at the lev e l of the bare NN interaction ( Amundsen and 0stgaard 



1985 



Baldo et al 



1992 



Elgar0y 



et a[], |l996a ; Kodcl et al , 1996; Takatsuka and Tamagaki , 1993| ). However, when comparing the results, one should 



note that the NN interactions used in these calculations are not phase-s hift equivalent, i.e., they do not predict exactly 



the same NN phase shifts. Furthermore, for the inte r action s used in ( Amundsen and 0stgaard , 1985| ; Baldo et al 



1992| ; |Elgar0y et al , 1996a; Takatsuka and Tamagaki, 1992) the predicted phase shifts do not agree accurately with 



modern phase shift analyses, and the fit of the NN data has typically x^/datum « 3. 

Fig. ^contains a comprehensive collection of our results for the ^P2-'^^2 pairing gaps with different potential models. 
We start with the top part of the figure that displays the results calculated with free single-particle energies. Differences 
between the results are therefore solely due to differences in the ^P2-^F2 matrix elements of the potentials. The plot 
shows results obtained with the old as well as with the modern potentials. The results (with the notable exception of 
the Argonne Vu interaction model) are in good agreement at densities below kp « 2.0 fm~^, but differ significantly 
at higher densities. This is in accordance with the fact that the diagonal matrix elements of the potentials are very 
similar below kp ~ 2.0 fm"^, corresponding to a laboratory energy for free NN scattering of i^iab ~ 350 MeV. This 
indicates that within this range the good fit of the potentials to scattering data below 350 MeV makes the ambiguities 
in the results for the energy gap quite small, although there is, in general, no unique relation between phase shifts 
and gaps. 

We would also like to calculate the gap at densities above kp ~ 2.0 fm^^. Then we need the various potentials at 
higher energies, outside of the range where they are fitted to scattering data. Thus there is no guarantee that the 
results will be independent of the model chosen, and in fact the figure shows that there are considerable differences 
between their predictions at high densities, following precisely the trend observed in the phase-shift predictions: the 



Argonne Vis is the most repulsive of t he modern potentials, followed by the CD-Bonn (Machleidt et al, 1996) and 
Nijmegen I and II (3toks et al, 1994). Most remarkable are the results obtained with Nijm-II: we find that the 



predicted gap continues to rise unrealistically even at kp « 3.5 fm~ , where the purely nuclconic description of 
matter surely breaks down. 

Since the potentials fail to reproduce the measured phase shifts beyond i^iab = 350 MeV, the predictions for the 
'^P2-^F2 energy gap in neutron matter cannot be trusted above kp « 2.0 fm"^. Therefore, the behavior of the ^P2-^F2 
energy gap at high densities should be considered as unknown, and cannot be obtained until potential models which 
fit the phase shifts in the inelastic region above Slab = 350 MeV are constructed. These potential models need the 
flexibility to include both the flat structure in the phase shifts above 600 MeV, due to the NN NA channel, as well 
as the rapid decrease to zero at Euh ~ HOO MeV. 

We proceed now to the middle part of Fig. |^, where the results for the energy gap u sing Brueckner-Hartr ec-Fock 



(BHF) single-particle energies are shown. For details on the BHF calculations, see, e.g., (Jeukenne et al 
this figure, two trends are apparent 
of the 3a 



1976). From 



First, the reduction of the i n-medium nucleon mass le a ds to a sizeable reduction 
F2 energy gap, as observed in earlier calculations ( Anmndsen and 0stgaar^ , 1985; Baldo et al, 1992 
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Elgar0y et al, 1996a; Takatsuka and Tamagaki , 1993). Secondly, the new NN interactions give again similar results 
at low densities, while beyond kp « 2.0 fm~ the gaps differ, as in the case with free single-particle energies. 

The single-particle energies at moderate densities obtained from the new potentials are rather similar, particularly 
in the important region near kp. This is illustrated by a plot. Fig. R of the neutron effective mass. 
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(36) 



as a function of density. Up to kp w 2.0 fm^ all results agree satisfactorially, but beyond that point the predictions 
diverge in the same manner as observed for the phase shift predictions. The differences in the BHF gaps at densities 
slightly above kp « 2.0 fm~^ are therefore mostly due to the differences in the ^P2-'^F2 waves of the potentials, but 
at higher densities the differences between the gap are enhanced by differences in the single-particle potentials. An 
extreme case is again the gap obtained with Nijm-II. It is caused by the very attractive matrix elements, amplified 
by the fact that the effective mass starts to increase at densities above kp « 2.5 fm^^ with this potential. 

Finally, in the lower panel of Fig. we illustrate the effect of different approximation schemes with an individual 
NN potential (CD-Bonn), namely we compare the energy gaps obtained with the free single-particle spectrum, the 
BHF spectrum, and an effective mass approximation, 



e(fc) = Uo 



fc2 



2m* ' 



(37) 



where m* is given in Eq. (|36D. In addition, also the gap in the uncoupled P2 channel, i.e., neglecting the tensor 
coupling, is shown. 

It becomes clear from the figure that the BHF spectrum forces a reduction of the gap by about a factor of 2-3. 
However, an effective mass aproximation should not be used when calculating the gap, because details of the single- 
particle spectrum around the Fermi momentum are important in order to obtain a correct value. The single-particle 
energies in the effective mass approximation are too steep near kp. We also emphasize that it is important to solve 
the coupled ^P2-^F2 gap equations. By eliminating the '^P2-^F2 and '^F2 channels, one obtains a ^P2 gap that is 
considerably lower than the ^P2-^F2 one. The reduction varies with the potential, due to different stren gths of the 
tens or force. For more detailed discussions of the importance of the tensor force, the reader is referr ed to (Amundsen 



and |0stgaardi |1985| ; |Elgar0y et al\ , |l996a| ; |Kodcl et al.[ |1998| , |200l| ; [lakatsuka and Tamagakij |1993D 



We end this subsection with a discuss ion of pairing for /3-stable matter of relevance for neutron star cooling, see for 

We will also omit a discussion on neutron pairing gaps in the ^Sq channel, 



example ( Pethick , 1992; Tsuruta 

since these appear at densities corresponding to the crust of the neutron star, see for example (|Barranco et 



1997) 



The gap in the crustal material is unlikely to have any significant effect on cooling processes ( Pethick and RavenhaH , 



1995), though it is expected to be important in the explanation of glitch phenomena. Therefore, the relevant pairing 
gaps for neutron star cooling should stem from the proton contaminant in the ^Sq channel, and superfiuid neutrons 
yielding energy gaps in the coupled ^P2-^F2 two-neutron channel. 

To obtain an effective interaction and pertinent single-particle energies at the BHF level, we can easily solve the 
BHF equations for different proton fractions. The conditions for f3 equilibrium require that 



where fii is the chemical potential of particle type i, and that charge is conserved 

TLr) — Tie 1 



(38) 



(39) 



where rii is the particle number density for particle i. If muons are present, the condition for charge conservation 
becomes 



and conservation of energy requires that 



He + Tin 



(40) 



(41) 



We assume that neutrinos escape freely from the neutron star. The proton and neutron chemical potentials are 
determined from the energy per baryon, calculated self-consistently in the MBHF approach. The electron chemical 
potential, and thereby the muon chemical potential, is then given by fie = fJ-n — fJ-p- The Fermi momentum of lepton 
type Z = e, /i is found from 



mf 



(42) 
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where m/ is the mass of lepton /, and we get the particle density using n; = kf /Stt"^. The proton fraction is then 
determined by the charge neutrahty condition (^). 

Since the relevant total baryonic densities for these types of pairing will be higher than the saturation density of 
nuclear matter, we will account for relativistic effects as well in the calculation of the pairing gaps. As an example, 
consid er the evaluation of the proton pairing gap using a Dirac-Brueckncr-Hartrcc-Fock approach, see (Elgar0y 
et a[], 1996a| ,|h|) for details. In Fig. ^ we plot as a function of the total baryonic density the pairing gap for protons 
in the ^Sq state, together with the results from a standard non- relativistic BCS approach. These results are all for 
matter in /3-equilibrium. In Fig. ^ we also plot the corresponding relativistic results for the neutron energy gap in the 
channel. For the ^Pq and the channels, the non-relativistic and the relativistic energy gaps vanish. 
As can be seen from Fig. ^ there are only small differences (except for higher densities) between the non-relativistic 
and relativistic proton gaps in the ^Sq wave. This is expected since the proton fractions (and their respective Fermi 
momenta) are rather small; however, for neutrons, the Fermi momenta are larger, and we would expect relativistic 
effects to be important. At Fermi momenta which correspond to the saturation point of nuclear matter, kp = 1.36 
fm^^, the lowest relativistic correction to the kinetic energy per particle is of the order of 2 MeV. At densities higher 
than the saturation point, relativistic effects should be even more important. Since we are dealing with very small 
proton fractions a Fermi momentum of kp = 1.36 fm~^ would correspond to a total baryonic density 0.09 fm~^. 
Thus, at larger densities, relativistic effects for neutrons should be important. This is also reflected in Fig. ^ for 
the pairing gap in the '^P2 channel. The maximum of the relativistic ^P2 gap is less than half the corresponding 



non-relativistic one and the density region over which it does not vanish is also much smaller; see (Elgar0y et 



1996b) for further details. 



This discussion can be summarized as follows. 



The 5*0 pro ton gap in /3-stable matter is < 1 MeV, and if polarization effects were taken into account (Schulze 



et ai, 1996| ), it could be further reduced by a factor of 2-3. 



• The ^P2 gap is also small, of the order of ^ 0.1 MeV in /3-stable matter. If relativistic effects are taken into 
account, it is almost vanishing. However, there is quite some uncertainty with the value for this pairing gap for 
densities above ^ 0.3 fm"'^ due to the fact that the NN interactions are not fitted for the corresponding lab 
energies. 



• Higher partial waves give essentially vanishing pairing gaps in /3-stable matter. 

Thus, the and ■^P2 partial waves are crucial for our understanding of superfluidity in neutron star matter. 

As an exotic aside, at densities greater than two-three times nuclear matter saturation density, model calculations 

iStoks and Lc^ , |2000| ; |Stoks and Rijkcn| , |1999| ; [Vidaiia 



based on baryon-baryon interactions ( p3aldo et ai , 1998a, 2000|; 



et Q^, |2000| ) or relativistic mean field calculations ( piendcnning , 2000| ) indicate that hyperons like S and A are 
likely to appear in neutron star matter. The size of the pairing gaps arising from these bar yons is , however, still 



an open problem, as it depends enti r ely o n the parametrization of the interaction models, see (Balberg and Barnea, 



1997 ; Schaab et ai , 1998 ; Takatsuka , 2002 ) for a critical discussion. Preliminary calc ulations of the pa i ring g ap for A 



hyperons using recent meson-exchange models for the hyperon-hyperon i nteraction (3toks and Rijkcn , 1999) indicate 
a vanishing gap, while S~-hyperon has a gap of the size of several MeVs ( Elgar0y and Schulze , 2001), At large baryon 
den sities f or wh ich perturbative QCD applies, pairing gaps for like quarks have been estimated to be a few MeV ( p] 
and Lovel, 1984), However, the pairing gaps of unlike q uarks {ud, us , and ds) have been suggested to be several tens 



to hundreds of MeV through non-perturbative studies (Alford et al, 1999), 

The cooling of a young (age < 10^ yr ) neutron star is mainly governed by neutrino emission processes and the specific 



heat (Page et ai, ^000 ; Schaab et ai, 1997, 1996), Due to the extremely high thermal conductivity of electrons, a 
neutron star becomes nearly isothermal within a time tu, ^ I — 100 years after its birth, depending upon the thickness 



of the crust (Pcthick and Ravcnhall, 1995), After this time, its thermal evolution is controlled by energy balance: 



dE, 



th 



dT 

dt 



(43) 



where Eth is the total thermal energy and Cy is the specific heat. and are the total luminosities of photons 
from the hot surface and neutrinos from the interior, respectively. Possible internal heating sources, due, for example, 
to the decay of the magnetic field or friction from differential rotation, are included in $. Cooling simulations are 
typically performed by solving the h eat transpo r t and hydrostatic equations including general relativistic effects, see 



for example the work of Page et al. (Page et al, 200C) 



The most powerful energy losses are expected to be given by the direct URCA mechanism 



(44) 
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However, in the outer cores of massive neutron stars and in the cores of not too massive neutron stars (M < 
1.3 — IAMq), the direct URCA process is allowed at densities where the momentum conservation kp<k^+ kp is 
fulfilled. This happens only at densities p several times the nuclear matter saturation density po = 0.16 fni^-^. 

Thus, for a long time the dominant proc esses for neutrino emission have been the modified URCA processes. See, 



for example, (Pethicli, 



1992 



Tsuruta, 1998) for a discussion, in which the two reactions 



P 



P 



(45) 



occur in equal numbers. These reactions are just the usual processes of neutron /3-decay and electron capture on 
protons of Eq. (^), with the addition of an extra bystander neutron. They produce ncutrino-antineutrino pairs, but 
leave the composition of matter constant on average. Eq. ( ^ ) is referred to as the neutron branch of the modified 
URCA process. Another branch is the proton branch 



n + p^p + p + e- 



P + P + e 



-p+Ve 



(46) 



Similarly, at higher densities, if muons are present, we may also have processes where the muon and the muon 
neutrinos (V^ and i'^) replace the electron and the electron neutrinos {Ve and i^e) in the above equations. In addition, 
one also has the possibility of neutrino-pair bremsstrahlung, processes with baryons more massive than the nucleon 
participating, such as isobars or hyperons or neutrino emission from more exotic states like pion and kaon condensates 
or quark matter. 

There are several cooling calculations including both superfluidity and many of the above processes, see for example 
(Page et ai, 2000; Schaab et ai, 1997, 1996). Bot h normal n e utron star matter and exotic states such as hyperons 



are included. The recent simulation of Page et al. ( Page et ai , 200C) seems to indicate that available observations of 
thermal emissions from pulsars can aid in constraining hyperon gaps. However, all these calculations suffer from the 
fact that the microscopic inputs, pairing gaps, composition of matter, emissivity rates, etc. are not computed at the 
same many-body theoretical level. This leaves a considerable uncertainty. 

These calculations deal however with the interior of a neutron star. The thickness of the crust and an eventual 
superfluid state in the crust may have important consequences for the surface temperature. The time needed for a 
temperature drop in the core to affect the surface temperature should depend on the thickness of the crust and on 
its thermal properties, such as the total specific heat, which is strongly influenced by the superfluid state of matter 
inside the crust. 

It has recently been proposed that the Coulomb-lattice structure of a neutron star crust may influence significantly 
the therm o dyna mical properties of the superfluid neutron gas (Broglia et ai, 1994). The authors of (Pcthick and 



Rav3nhal], 1995| ) have proposed that in the crust of a neutron star non-spherical nuclear shapes could be present 



at densities ranging from p — 1.0 x 10^** gcm"'^ to p = 1.5 x lO^'' gem"'', a density region which represents about 
20% of the wh ole crust. The saturation den sity of nuclear matter is po = 2.8 x 10^^ gcm~'^. These unusual shapes 
are supposed ( Pcthick and Ravenhall , 1995 ) to be disposed in a Coulomb lattice embedded in an almost uniform 
background of relativistic electrons. According to the fact that the neutron drip point is supposed to occur at lower 
density (p ~ 4.3 x 10^^ gcm"'^), and considering the characteristics of the nuclear force in this density range, we 
expect these unusual nuclear shapes to be surrounded by a gas of superfluid neutrons. 



To model the influence on the heat conduction due to pairing in the crust, Broglia et al (Broglia et al, 1994) studied 



various nuclear shapes for nuclei immersed in a neutron fluid using phenomenological interactions and employing a 
local-density approach. They found an enhacement of the fermionic specific heat due to these shapes compared to 
uniform neutron matter. These results seem to indicate that the inner part of the crust may play a more relevant 
role on the heat diffusion time through the crust. Calculations with realistic nucleon- nucleon interaction were later 
repeated by Elgar0y et al (Elgar0y et ai, 1996d), with qualitatively similar results. 



2. Proton-neutron pairing in symmetric nuclear matter 



The calculation of the ^5*0 gap in symmetric nuclear matter is closely related to the one for neutron matter. Even 
with modern charge-dependent interactions, the resulting pairing gaps for this partial wave are fairly similar, see for 
example ( Elgar0y and Hjorth- Jensen] , 1998). 

The size of the neutron-proton (np) '^Si-^Di energy gap in symmetric or asymmetric nuclear matter has, however, 
been a much debated issue since the first calculations of this quantity appeared. While solutions of the BCS equations 
with bare nucleon- nucleon (NN) forc e s give a large energy gap of several MeVs at the saturation d e nsity kp = 1.36 fm ~^ 



(p ^ 0.17 fm_^ (|Alm et al. 



1990 



Baldo et al 



1995 



Garrido et al. 



Lonbardc, 2000; Takatsuka and Tamagaki , 1993; Vonderfecht et al, 1991D , there is little empirical evidence from finite 
nuclei for such strong np pairing correlations, except possibly for isospin T = and N = Z, see also the discussion in 



2001 



Scdrakian et al, 1997: Sedrakian and 
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Sec. Ill and the recent work of Jenkins et al (Jenkins et al, 2002). One possible resolution of this problem lies in the 



fact that all these calculations have neglected contributions from the induced interaction. Fluctuations in the isospin 
and the spin-isospin channel will probably make the pairing interaction more repulsive, leading to a substantially 
lower energy gap. One often-neglected aspect is that all non-relativistic calculations of the nuclear matter equation 
of state (EOS) with two-body NN forces fitted to scattering data fail to reproduce the empirical saturation point, 
seemingly regardless of the sophistication of the many-body scheme employed. For example, a BHF calculation of 
the EOS with recent parametrizations of the NN interaction would typically give saturation at kp = 1.6-1.8 fm^^. 
In a non-relativistic approach, it seems necessary to invoke three-body forces to obtain saturation at the empirical 
equilibrium density, see for example (Akmal et al, 199S). This leads one to be cautious when talking about pairing at 



the empirical nuclear matter saturation density when the energy gap is calculated within a pure two-body force model, 
as this density will be below the calculated saturation density for this two-body force, and thus one is ca lculating 
the gap at a density where the system is theoretically unstable. One even runs the risk, as pointed out in (Jackson, 
1983| ) , that the compressibility is negative at the empirical saturation density, which means that the system is unstable 



against collapse into a non-homogeneous phase. A three-body force need not have dramatic consequences for pairing, 
which, after all, is a two-body phenomenon, but still it would be of interest to know what the ^Si-^Di gap is in a model 
in which the saturation properties of nuclear matter are reproduced. If one abandons a non-relativistic description, 
the empirical saturation point can be obta ined within the Dirac-Brueckner- Hartree-Fock (DBHF) approach, as first 
pointed out by Brockmann and Machleidt (Brockmann and Machleidt, 1990). This might be fortuitous, since, among 
other things, important many-body effects are neglected in the DBHF approach. Nevertheless, it is interesting to 
investigate '^Si-^Di pairing in this model and compare our results with a corresponding non-relativistic calculation. 
Furthermore, several groups have recently developed relativistic formulations of pairing in nuclear matter (Guimaraes 
et a[], \199t ; Kucharek and Ring, 1991;|Matera et al. , 1997 ;_ Serra et al , p002 ) and have applied them to ^5*0 pairing. 
The models are of the Walecka-type (Berot and Walecka, 1986| ) in the sense that meson masses and coupling constants 
are fitted so that the mean-field EOS of nuclear matter meets the empirical data. In this way, however, the relation of 



the models to free-space NN sca t tering becomes som ewhat unclear. An interesting result found in (Guimaraes et al 



1996: 



Kucharek and Ring , 



1991; 



Matera et al\ |l997| ) is that the ^Sq energy gap vanishes at densities slightly below 



the empirical saturation density. This is in contrast with non-rela tivistic calcu l ation s which generally give a relatively 



small, but non- vanishing ^Sp gap at this density, see for instance ( Baldo et al , 199C; Chen et al , 1982; Elgar0y et al 



1996c; Kucharek et al, 1989) 



In Fig^ |10| we show the EOS obtained in our non-relativistic and relativistic calculations. The non-relativistic 
one fails to meet the empirical data, while the relativistic calculation very nearly succeeds. In these calculations, 
we employed the non-relativistic and relativistic one-boson exchange models from the Bonn A interaction defined 
in (Machleidt, 1989). A standard BHF calculation was done in the non-relativistic case, whereas in the relativistic 
case we incorporate minimal relativity in the gap equation, thus using DBHF single-particle energies in the energy 
denominators and modifying the free NN interaction by a factor fh? / EkE^' (Elgar0y et al, 1996b). The resulting 
pairing gaps are shown in Fig. |ll|. For the non-relativistic calculation, we see a large energy gap at the empirical 
saturation density around 6 MeV at kp — 1.36 fm~^, in agreement with earlier non-relativistic calculations (Aim et al 



199C; Baldo et al, 1995; Takatsuka and Tamagaki, 1993; Vondcrfecht et al, 1991). In the relativistic calculation, we 



find that the gap is vanishingly small at this density. 

Since non-relativistic calculations with two-body interactions will, in general, give a saturation density that is too 
high (an example is shown in Fig. [l0| ) , this implies that in a non-relativistic approach we are actually calculating the 
gap at a density below the theoretical saturation density, and one may question the physical relevance of a large gap at 
a density where the system is theoretically unstable. If one considers the gap at the calculated saturation density for a 
non-relativistic approach with a two-body force only, it is in fact close to zero. In the DBHF calculation, we come very 
close to reproducing the empirical saturation density and binding energy, and when this is used as a starting point for 
a BCS calculation, we find that the gap vanishes, both at the empirical and the calculated saturation density. That the 
DBHF calculation meets the empirical points is perhaps fortuitous, as important many-body diagrams are neglected 
and only medium modifications of the nucleon mass are accounted for. An increased repulsion in the non-relativistic 
may thus reduce the gap dramatically. 

We end this section with a comment on the interesting possibility of a transition from BCS pairing to a Bose- 
Einstein condensantion in asymmetric nuclear matter at low densities. For the singlet ^Sq partial wave we do not 
expect to see a transition, essentially because the coherence length is much bigger than the interparticle spacing. The 
inclusi on of medium effects such as screening terms are expected to further reduce the pairing gap, see (De Blasio 
et a[], 1997 ) and thereby enhance the coherence length. However, this does not imply that such a transition is not 
possibl e in nuclear matter or asymm etric nuclear matter as present in a neutron s t ar. A recent analysis by Lombardo 



et al (Lombardo and Schuck, 2001), see also the work of Baldo et al. (Baldo et al, 1995), of triplet ^Si pairing in low- 



density symmetric and asymmetric nuclear matter, indicates that such a transition is indeed possible. As the system 
is diluted, the BCS state with large overlapping Cooper pairs evolves smoothly into a Bose-Einstein condensation 



19 



of tightly bound deuterons, or neutron-proton pairs. A neutron excess in this low-density regime does not affect 
these deuterons due to the large spatial separation of the deuterons and neutrons. Even at large asymmetries, these 
deuterons are only weakly affected. This effect can have interesting consequences for the understanding of e.g., exotic 
nuclei and asymmetric and expanding nuclear matter in heavy-ion collisions. 



E. Conclusions and open problems beyond BCS 



We have seen that pairing in neutron star matter is essentially determined by singlet pairing in the ^Sq channel 
and triplet pairing in the channel. These two partial waves exhibit a contribution to the NN interaction which 
is attractive for a large range of densities. These partial waves are also crucial for our understanding of pairing 
correlations in finite nuclei. Whether it is possible to have a strong neutron-proton pairing gap for symmetric matter 
in the '^Si channel is still an open question. Relativistic calculations indicate a vanishing gap at nuclear matter 
saturation density. The results we have discussed have all been within the frame of a simple many-body approach; 
however, the analyses that have been pe r formed are not co ntingent upon these simplifications. Combined with, e.g. 



the separation analysis of ( Kodel et a/. , 1996 , 1998 , 2001), we believe the calculation procedures will retain their 



validity when more complicated many-body terms are inserted. 

A complete and realistic treatment of pairing in a given strongly coupled Fermi system such as neutron matter 
demands ab initio calculation of both the single-particle energies and the interaction in the medium. The dependence 
of, e.g., pairing upon various approaches to the single-particle energies is a clear signal of the need for a consistent 
many-body scheme, see for example Fig. 0. Whether we employ a density-dependent effective mass approach as in 
Eq. (36) or a standard effective mass approach as in Eq. (p7|), the results are dif ferent contributions to the pairing 



gap. Recently, Lombardo et al. (Lombardo et al, 2001b; Lombardo and Schulze, 2001| ) reexamined the role played 
by ground-state correlations in the self-energy. Solving the Gorkov equations, see Eq. (|^), they found a substantial 
suppression of the ^5*0 pairing due to changes in the quasiparticle strength around the Fermi surface. Their results 
are shown in Fig. |l2| for a set of different /c^-values. 

This figure shows that self-energy effects are an important ingredient in our understanding of the pairing gap in 
infinite matter. 

A correct treatment of the self-energy entails a self-consistent scheme where the renormalization of the interaction is 
done at an equal footing. Of special interest for the pairing interaction are polarization corrections. At low densities we 
may expect that the dominant polarization term stems from a second-order perturbative correction with particle-hole 
intermediate states, as depicted in Fig. |l^. 

For contributions around the Fermi surface, one can evaluate diagram (a) analytically and obtain a result in terms 
of the Ferm i momentum and th e scattering length. As shown by Heiselberg et al. and Schulze et al. (Heiselberg 
et a I, |2000| ; ^chulze et al| , ^OOlj ), even the low-density expression of Eq. (E9[) is reduced by a factor of « 2.2 when 
polarization terms are included. 

To go beyond diagram (a) and simple low-density approximations requires considerable efforts and has not been 
accomplished yet. This means that there is still a large uncertainty regarding the value of the pairing gap in infinite 
matter. There are few calculations of the pairing gap from the point of view of an ab initio approach 



On e such scheme is the on e favored by Clark and co-workers, based on correlated-basis (or CBF) theory ( Bishop , 



1991; Chen et al, 1986| , 1983). Within the CBF scheme, the following approach to the quantitative physics of pairing 



in extended nucleonic systems has been undertaken: 



(a) Dressing of the pairing interaction by Jastrow correlations within CBF theory ( [Krotscheck and Clark] , |1980 



Krotscheck et 



1981) 



(b) Dressing of the pairing interaction by dynamical collective effects within CBF theory ( [Chen et al , 198(;, 1983 



Krotscheck et al 
ations, etc.) 



1981) (including polarization effects arising from exchange of density and spin-density fluctu- 



(c) Consistent renormalization of single-particle energies by short- and long-range correlations within CBF theory 
(cf. ( [Krotscheck and Clark) , |1983| )) 



This approach has already been explored in the ^So neutron pairing problem ( phen et al , 1986, 1983) ), although 



the assumed Jastrow correlations have not been optimized and only a second-order CBF perturbation treatment is 
available for step (b). Application of this scheme t o •^P2 ~'^^2 pairing in neutron-star matter is still an unexplored 
topic. Alternatively, coupled-cluster (CC) (Bishop, 1991) or Fermi-hypernetted chain inspired approaches could be 



used (Fabrocini et al, 1998). Another app r oach, followed by Wambach, Ainsworth and Pines (Ainsworth et al 



198£, 1993) and Schulze et al (Schulze et al, 1996|) d eparts from the Landau theory inspired many 



3ody approach to 





Dickhoff et al. 


1983, 


1981 


; Dickhoff and 
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Miither, 1987; Jackson et al, 1982a| ). This microscopic derivation of the effective interaction starts from the following 



physical idea: the particle-hole (p-h) interaction can be considered as made of a direct component containing the 
short-range correlations and an induced component due to the exchange of the collective excita tions of the medium. 



Finally, another alt ernative is to solve the full set of the Parquet equations, as discussed in ( Hjorth- Jensen , 



2002 



Jackson et al, 1982b). This self-consistent scheme entails the summation to all orders of all two-body diagrams with 



particle-particle and hole-hole (ladder diagrams) and particle-hole (polarization and screening diagrams) inter mediate 
state s, accompanied with the solution of Dyson's equation for the single-particle propagator. Recently, Bozek (Bozek, 
2002) has studied the generalized ladder diagram resummation in the superfluid phase of nuclear matter. This is the 
first step towards the solution of the Parquet diagrams. 

We conclude by summarizing this section through Fig. |lj. This figure exhibits the influence of various approaches 
which include screening corrections to the pairing gap. The curve in the background is given by the calculation 
with free single-particle energies and the bare nucleon-nucleon interaction. These calculations are similar, except for 
the potential model employed, to those discussed in, e.g., Fig. ^. This means that the calculations of Subsec. II. C, 
with only experimental inputs, phase shifts and scattering length, yield an upper limit for the pairing gap. How 
such renomalizations will affect the '^P2 gap is an entirely open issue. This gap is crucial since it extends to large 
densities and can reasonably be expected to occur at the centers of neutron stars. Unfortunately, one cannot constrain 
at p resen t the size of the pairing gap from data on thermal emission from neutron stars, see also the discussion in 
Sec. fflDl. 



III. PAIRING CORRELATIONS IN FINITE NUCLEI 
A. Introduction to the nuclear shell model 

Our tool for analyzing pairing correlations in finite nuclei is the nuclear shell model, with appropriately defined 
model spaces and effective interactions. In this section we extract information on pairing correlations through large- 
scale shell-model calculations of several nuclear systems, from nuclei in the sd-shell to heavy tin isotopes. 

We define the nuclear shell model by a set of spin-orbit coupled single-particle states with quantum numbers 
Ijm denoting the orbital angular momentum (l) and the total angular momenta (j) and its z-component, m. In 
a rotationally invariant basis, the one-body states have energy sij that are independent of m. The single-particle 
states and energies may be different for neutrons and protons, in which case it is convenient to include also the 
isospin component tz = ±1/2 in the state description. We will use the label a for the set of quantum numbers 
Ijm or Ijmtz, as appropriate. These orbits define the valence P-space, or model space for the shell model, while 
remaining single-particle orbits define the so-called excluded space, or Q-space. We can express these spaces through 
the operators 

n oo 

P = ^|V,)(V'd, Q= E (47) 

i—1 i—n-\-l 

where n defines the dimension of the model space while the wave functions ipi could represent a many-body Slater 
determinant built on the chosen single-particle basis. As an example, if we consider the chain of tin isotopes from 
^•^"Sn to ^■^^Sn, the neutron single-particle orbits 2si/2, 1^5/2: 1^^3/2i 0^7/2, and 0/iii/2 could define an eventual model 
space. We could then choose ^^^Sn as a closed-shell core. Neutron holes from ^^^Sn to ^'^"Sn define then the valence- 
space or model-space degrees of freedom. We could, however, have chosen ^°"Sn as a closed-shell core. In this case, 
neutron particles from ^°^Sn to ^'^^Sn define the model space. 

The shell-model Hamiltonian H is thus built upon such a single-particle basis. The shell- model problem requires 
normally the solution of a real and symmetric n x n matrix eigenvalue equation 

H\^'k)^Ek\^k), (48) 

with k = 1, . . . , n, where the size of this matrix is defined by the actual shell- model space. The dimensionality n of 
the eigenvalue matrix H is increasing with an increasing number of valence particles or holes. As an example, for 
^^^Sn with the above mentioned single-particle basis, the dimensionality of the Hamiltonian matrix is of the order of 
n ~ 10^. For nuclei in the rare-earth region, this dimensionality can be of the order of n ~ 10^^ — 10^''. 

The shell-model Hamiltonian can be written in the form H = Hi + H2 + H3 + ■ ■ ■ where Hi is a one-bod y t erm 
typically represented by experimental single-particle energies, see Eq. (20). The two-body term, see Eq. (po|), is 



given in terms of the uncoupled matrix elements V of the two-body interaction. These matrix elements must obey 
rotational invariance, parity conservation, and (when implemented) isospin invarience. To make explicit the rotational 
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and isospin invariance, we rewrite the two-body Hamiltonian as 



where the pair operator is 



^JT-MT, ("^) = E Uamajpmfi \ J M) ( ^tai^^g | TT^ j a] 



2 "2 



(50) 



In these expressions [JM) are the coupled angular momentum quantum numbers and (TT^) are the coupled isospin 
quantum numbers. The coupled two-body matrix elements Vjt define the valence particle interactions within the 
given shell- model space. They are matrix elements of a scalar potential V{'ri, r2) and are defined as 



(51) 



and arc independent of M and Tj. The antisymmetrized matrix elements are Vj^{aP,jS) and are then given by 

F/^(a/3,7<^) = [(l + Scp) {l + S,s)r'^' [VMaP,j6) - (-l)'^+^"+-''^+^- VjT(/?a, 7'^)] ■ (52) 
We remark here that three-body or higher-body terms such as are normally not included in a shell-model effective 



Miither et al. 



interaction, although shell- model analyses with three-body interactions have been made in (Engcland et al\, 2002 
|1985| ) 



In the following subsections, we discuss how to extract information about pairing correlat ions w ithin the framework 
of large-scale shell- model and shell-model Monte Carlo (SMMC) calculations. In Subsec. III.B, we discuss selected 
features of the tin isotopes such as the near constancy of the energy difference between the first excited state with 
J = 2 and the ground state with J = for the whole chain of even isotopes from ^^^Sn to ^'^"Sn. These are nuclei 
whose excited states are well reproduced by the neutron model space mentioned above. We relate this near constancy 
to strong pairing correlations and the same partial waves which contribute to superfiuiditity in neutron stars, namely 
the ^Sq and components of the nucleon-nucleon interaction. The ^Sq component is generally the dominating 
partial wave, a well-known fact in nuclear physics. We show also that a truncation scheme like generalized seniority 
(Talmi, 1993| ) is a viable first approximation to large-scale shell- model calculations. 



In Subsec. pLC 



we discuss isoscalar and isovector pairing correlations, whereas proton- neutron pairing and Wigner 
energy are discussed in Subsec. [II. D. Various thermal properties are discuss ed in the rema i ning subsections. These 
results are obtained through large-scale SMMC calculations, see for example (Koonin et al., 1997). 



B. Tin isotopes, seniority, and the nucleon-nucleon interaction 

Nuclei far from the line of /3-stability are at present in focus of the nuclear structure physics community. Considerable 
attention is being devoted to the experimental and theoretical study of nuclei near ^^''Sn from studies of the chain of 
Sn isotopes up to ^^^Sn to, e.g., nuclei near the proton drip line hke ^^^^^^^Sb. 

Our scheme to obtain an effective two-body interaction for shell-model studies starts with a free nucleon-nucleon 
interaction V , which is appropriate for nuclear physics at low and intermediate energies. Here we employ the charge- 
dependent version of the Bonn potential models, see (Machleidt, 2001) and the discussion in Sec. |l| The next step in 
our many-body scheme is to handle the fact that the repulsive core of the nucleon-nucleon potential V is unsuitable for 
perturbative approaches. This problem is overcome by introducing the reaction matrix G, which in a diagrammatic 
language represents the sum over all ladder types of diagrams. This sum is meant to renormalize the repulsive short- 
range part of the interaction. The physical interpretation is that the particles must interact with each other an infinite 
nu mber of times in order t o pro duce a finite interaction. We calculate G using the double-partitioning scheme discussed 



in (Hjorth- Jensen et al., 1995). Since the G-matrix represents just the summation to all orders of particle-particle 
ladder diagrams, there are obviously other terms which need to be included in an effective interaction. Long-range 
effects represented by core-polarization terms are also needed. In order to achieve this, the G-matrix elements are 
renormalized by the so-called Q-box method. The Q-box is made up of non-fol ded diagrams which are irre ducible and 



valence-linked. Here we include all non-folded diagrams to third order in G ( Hjorth- Jensen et al. , 1995 ). Based on 



the Q-box, we compute an effective interaction H in terms of the Q-box using the folded-diagram expansion method 
(see for example (Hjorth- Jensen et al., 1995) for further details). 
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The effective two-particle interaction is then used in larg e-scale sheh-model ca lculations. For the shell-model 
calculation, we employ the Oslo m-scheme shell-model code ( Engeland et al , 2002 ), which is based on the Lanczos 
alg orithm, an iterative met hod which gives the solution of the lowest eigenstates. The technique is described in detail 



in ( Whitehead et al. , 1977 ). The shell-model space consists of the orbits 2si/2, 1^5/21 l<^3/2i O57/2 ^nd O/i 



11/2- 



Of interest in this study is the fact that the chain of even tin isotopes from ^"^Sn to -'^'^"Sn exhibits a near constancy of 
the 2i — excitation energy, a constancy which can be related to strong pairing correlations and the near degeneracy 
in energy of the relevant single-particle orbits. As an example, we show the experimental 2^ — 0^ excitation energy 
from ^^^Sn to ^^''Sn in Table ||. Our aim is to see whether partial waves which play a crucial role in superfluidity of 
neutron star matter, viz. ^Sq and '^^2, are equally important in reproducing the near-constant spacing in the chain 
of even tin isotopes shown in Table |. 

In order to test whether the ^5*0 and partial waves are equally important in reproducing the near constant 
spacing in the chain of even tin isotopes as they are for the superfluid properties of infinite neutron star matter (recall 
the discussion of Sec. ||), we study four different approximations to the shell- model effective interaction, viz., 

1. Our best approach to the effective interaction, VcS_, contains all one-body and two-body diagrams through third 
order in the G- matrix, as discussed above, see also (Holt et al., 199S). 



2. The effective interaction is given by the G-matrix only and inludcs all partial waves up to Z = 10. 

3. We define an effective interaction based on a G-matrix which now includes only the partial wave. 

4. Finally, we use an effective interaction based on a G-matrix which does not contain the ^Sq and partial 
waves, but all other waves up to Z = 10. 



In all four cases the same NN interaction is used, viz., the CD-Bonn interaction described in ( Machleidf , 2001). Table 
I lists the results. 

We note from this table that the three first cases nearly produce a constant 2^ — Oj^ excitation energy, with our 
most optimal effective interaction Vcs being closest to the experimental data. The bare G-matrix interaction, with 
no folded diagrams as well, results in a slightly more compressed spacing. This is mainly due to the omission of the 
core-polarization diagrams which typically render the J = matrix elements more attractive. Such diagrams are 
included in Vcs- Including only the ^5*0 partial wave in the construction of the G-matrix (case 3) yields, in turn, a 
somewhat larger spacing. This can again be understood from the fact that a G-matrix constructed with this partial 
wave only does not receive contri butions from any entirely repulsive partial wave. It should be noted that our optimal 



interaction, as demonstrated in (Holt et al, 1998), shows a rather good reproduction of the experimental spectra for 



both even and odd nuclei. Although the approximations made in cases 2 and 3 produce an almost constant 2^ — 0^ 
excitation energy, they reproduce poorly the properties of odd nuclei and other excited states in the even Sn isotopes. 

However, the fact that the first three approximations result in a such a good reproduction of the 2^ — 0^ spacing 
may hint to the fact that the partial wave is of paramount importance. If we now turn attention to case 4, i.e., 
we omit the ^5*0 and partial waves in the construction of the G-matrix, the results presented in Table | exhibit 
a spectroscopic catastrophe^. We also do not list eigenstates with other quantum numbers. For ^^^Sn, the ground 
state is no longer a 0+ state; rather it carries J = 4+ while for ^^^Sn the ground state has 6+. The first 0^ state 
for this nucleus is given at an excitation energy of 0.1 MeV with respect to the 6+ ground state. The general picture 
for other eigenstates is that of an extremely poor agreement with data. Since the agreement is so poor, even the 
qualitative reproduction of the 2'^ — 0^ spacing, we defer from performing time-consuming shell-model calculations 

116418, 120422gj^ 

Since pairing is so prominent in such systems, we present a comparison of the SM with the generalized seniority 
model ( [Talmi , 1993). The generalized seniority scheme is an extension of the seniority scheme, i.e., from involving only 
one single j-orbital, the model is generalized to involve a group of j-orbitals within a major shell. The generalized 
seniority scheme is a more simple model than the shell model since a rather limited number of configurations with a 
strictly defined structure are included, thus allowing a more direct physical interpretation. States with seniority w = 
are by definition states where all particles are coupled in pairs. Seniority w = 2 states have one pair broken, seniority 
w = 4 states have two pairs broken, etc. The generalized seniority scheme is suitable for describing semi-magic nuclei 
where pairing plays an important role. The pairing picture and the generalized seniority scheme have been important 
for the description and understanding of the tin isotopes. A typical feature of the seniority scheme is that the spacing 



^ We limit the discussion to even isotopes from ^^''Sn to ^''"Sn, since a qualitatively similar picture is obtained from ^"^Sn to ^^®Sn. 
* Although we have singled out these two partials waves, due to their connection to infinite matter, it is essentially the ^ So wave which 
is responsible for the behavior seen in Table hi. 
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of energy levels is independent of the number of valence particles. For the tin isotopes, not only the spacing between 
the ground state and the 2^ state, but also the spacing beween the ground state and the 4^ and 6^ states is fairly 
consta nt thro u ghout the whole sequence of isotopes. In fundamental works on generalized seniority by, for instance, 
Talmi ( Talmi , 1993 ), the tin isotopes have been used as one of the major test cas e s. It is also worth mentioning 
t he classical work on pa iring by Kisslinger and Sorensen ( Kisslingcr and Sorcnsen , see also the analyses of 

( ^andulescu et al , 1997 ) and the review article of Bes and Sorensen ( Bes and Sorensenfigsg ). 

If we, by closer investigation and comparison of the SM wave function and the seniority states, find that the most 
important components are accounted for by the seniority scheme, we can benefit from this and reduce the SM basis. 
This would be particularly useful when we want to do calculations on systems with a large number of valence particles. 

The operator for creating a generalized seniority (w = 0) pair is 



5t 



^y2jTT' 



(53) 



m>0 



where fcj^ is the creation operator for holes. The generalized senitority {v = 2) operator for creating a broken pair is 
given by 



Ed 

j<j' 



''^^Pj,r{jmj'fn'\JM)b]^b],^, 



(54) 



The coefficients aj and f3jj' are obtained from the ^""^Sn ground state and the excited states, respectively. 

We calculate the squared overlaps between the constructed generalized seniority states and our shell-model states 



{v = 0) |(-4Sn(SM);0+|(5t)f |0)|2, 

iv^2) |(^Sn(SM);J,|7^t^,(5t)f-i| 



(55) 



The vacuum state |0) is the ^'^^Sn-core and n is the number of valence particles. These quantities tell to what extent 
the shell-model states satisfy the pairing picture, or in other words, how well is generalized seniority conserved as a 
quantum number. 

The squared overlaps are tabulated in Table and vary generally from 0.95 to 0.75. As the number of valence 
particles increases, the squared overlaps gradually decrease. The overlaps involving the 4+ states show a fragmentation. 
In ^2*Sn, the A'^ (SM) state is mainly a seniority v = 2 state. As approaching the middle of the shell, the next state, 
4^, takes more and more over the structure of a seniority v = 2 state. The fragmentation of seniority over these two 
states can be understood from the fact that they are rather close in energy and therefore may have mixed structure. 

In summary, these studies show clearly the prominence of pairing correlations in nuclear systems with identical 
particles as effective degrees of freedom. There is a clear link between superfluidity in infinite neutron star matter 
and spectra of finite nuclei such as the chain of tin isotopes. This link is provided especially by the ^5*0 partial wave of 
the nucleon-nucleon interaction. Excluding this component from an effective interaction yields spectroscopy in poor 
agreement with experimental data. Although the ^5*0 partial wave plays an important role, other many-body effects 
arising e.g., from low- lying collective surface vibrations among nucleons can have important effects on properties of 



nuclei, as demonstrated recently by Broglia et al (Barranco et al, 1999; Giovanardi et al., 2002). In order to interpret 
the results of Table | one needs to analyze the core-polarizations diagrams which are used to compute the effective 
interaction in terms of the various partial waves. 

Generalized seniority provides an explicit measure of the degree of pairing correlations in the wave functions. 
Furthermore, generalized seniority can serve as a useful starting point for large-scale shell-model calculations and is 
one among several ways of extracting information about pairing correlations. In the next subsections, we present 
further approaches. 



C. Isoscalar and isovector pair correlations 

Numerous phenomenological descriptions of nuclear collective motion describe the nuclear ground state and its 
low-lying excitations in terms of bosons. In one such model, the Interacting Boson Model (IBM), L = (S) a nd L = 2 
(D) bosons are identified with nucleon pairs having the same quantum numbers (lachello and Arima, 1988| ), and the 
ground state can be viewed as a condensate of such pairs. Shell- model studies of the pair structure of the ground 
state and its variation with the number of valence nucleons can therefore shed light on the validity and microscopic 
foundations of these boson approaches. 

Generally speaking, nucleon-nucleon pairing may be considered in several classes. A nucleon has a spin j — 1/2, jz — 
±1/2 and an isospin t = 1/2, t^ = ±1/2. Two protons (neutrons) thus are allowed to become paired to total J, T = 0, 1 
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and = —1 (T^ = 1). We shall call this isovector pairing. Isoscalar pairing delineates proton-neutron pairing for 
which J, T = 1, and = 0. 

While we concentrate here on shell model results, we do wish to point out to the reader several recent developments in 
other model studies of nuclcon-nucleon pairing. Interesting studies of nuclcon-nucl con pairing have b e en un dertaken 



in several models, including pscudo-S'[/(4) symmetry studies for p/-shcll nuclei (Van Isacker et al, 1999). These 



studies indicated that pseudo-iSt/ (4) is a reasonable starting point for the description of systems within the p/-shell 
larger than ^^Ni. It is also the starting point for generating collective pairs within th e framework o f the Interacting 
Boson Model that incorporates T = and T — 1 bosons and a bosonic SU{4) algebra ( Elliott , 1958 ). This symmetry 
dictates that pairing strengths are the same in the both the T — and T = 1 channels. Extensiv e studies of pairin g 
in the framework of Hartree-Fock-Bogoliubov theory have also been undertaken (see, for example (Goodman, 2000)). 
Recent work in this direction indicates that T — and T = 1 pairing superfluids may develop near the mid-point of 
isotope chains (i.e., near N = Z nuclei). 

In the framework of the shell model, it appears sufficient for many purposes to study the ECS pair structure in the 
ground state. The BCS pair operator for protons can be defined as 



At 



Et t 
PjmP 



(56) 



where the sum is over all orbitals with m > and — (— )■' "^p]m time-reversed operator. Thus, the 

observable A^A and its analog for neutrons are measures of the numbers of J = 0,T = 1 pairs in the ground state. 
For an uncorrelated Fermi gas, we have simply 



(At A) = y: 



2{2j + 1) 



(57) 



where the nj = {PjmPjm) are the occupation numbers, so that any excess of (At A) in our SMMC calculations over 
the Fermi-gas value indicates pairing correlations in the ground state. 

In this analysis, we move from tin isotopes to nuclei which can be described by the single-particle orbits of the 
p/-shell, Ijp.s/ T, IjPi /s^, O/s/9 , and O /7/9. As effective interaction, we employ the phenomenological interaction of Brown 
and Richter ( Richter et al. , 1991). Fig. [ij shows the SMMC expectation values of the proton and neutron BCS-like 
pairs, obtained after subtraction of the Fermi gas value (Eq. for three chains of isotopes. As expected, these excess 
pair correlations are quite strong and reflect the well-known coherence in the ground states of even-even nuclei. Note 
that the proton BCS-like pairing fields are not constant within an isotope chain, showing that there are important 
proton-neutron correlations present in the ground state. The shell closure at A'^ = 28 is manifest in the neutron 
BCS-like pairing. As is demonstrated in Fig. the proton and neutron occupation numbers show a much smoother 
behavior with increasing A. 

It should be noted that the BCS form (Eq. p6[ ) in which all time-reversed pairs have equal amplitude is not necessarily 
the optimal one a nd allows only the study of S-pair structure . To explore the pair content of the ground state in a 
more general way ( Alhassid et al. , 1996 ; Langanke et al. , 1996 ), we define proton pair creation operators 



(58) 



These operators are boson-like in the sense that 



[A]j^,{ja]b), Aj^iiajb)] = 1 + 0{h/2j + 1) ; 



(59) 



i.e., they satisfy the expected commutation relations in the limit of an empty shell. We may also construct from these 
operators a pair matrix 



1 



, = — , 



^^^/M Ua,jb)AjM (jc, id)) 



(60) 



We construct bosons B'^^j^ 



as 



^^a\{jajb)Al^{jajb) , 

jajb 



(61) 



25 



where a = 1, 2, . . . labels the particular boson and the "wave function" ip satisfies 



Yi'*ajUajb)i'pj{jajb) 



0af3 



(62) 



Jajb 



(Note that V' is independent of /i by rotational invariance.) 

To find and Uaj = '^fj{Bl^j^Bajfj.), the number of bosons of type a, and multipolarity J, we compute the 
quantity '^^{A''j^{jajb)^jfj.{jcjd)), which can be thought of as an hcrmitian matrix M'^^, in the space of orbital pairs 
Uajb)', its non-negative eigenvalues define the n^j (we order them so that nu > n2j > ■ ■ ■), while the normalized 
eigenvectors are the ^aj{jajb)- The index a distinguishes the various possible bosons. For example, in the complete 
p/-shell the square matrix M has dimension Nj = 4 for J = 0, Nj — 10 for J — 1, Nj = 13 for J = 2, 3. 

The presence of a pair condensate in a correlated ground state will be signaled by the largest eigenvalue for a given 
J, nij, being much greater than any of the others; Tpu will then be the condensate wavefunction. In Fig. |l^ we show 
the pair matrix eigenvalues n„j for the three isovector J — pairing channels as calculated for the iron isotopes 
54-58pg^ We compare the SMMC results with those of an uncorrelated Fermi gas, where we can compute {A^ A) using 
the factorization 



npUaiSfj^SaS - 5l3s5a-i) , (63) 

a'pCLp) are the occupation numbers. Additionally, Fig. |l^ shows the diagonal matrix elements of the 
As expected, the protons occupy mainly f'j/2 orbitals in these nuclei. Correspondingly, the {A^ A) 



where the = 
pair matrix Af, 

expectation value is large for this orbital and small otherwise. For neutrons, the pair matrix is also largest for the 
orbital. The excess neutrons in ^^-^^Fe occupy the P3/2 orbital, signaled by a strong increase of the corresponding pair 
matrix element M22 in comparison to its value for ^^Fe. Upon closer inspection, we find that the proton pair matrix 
elements are not constant within the isotope chain. This behavior is mainly caused by the isoscalar proton-neutron, 
pairing. The dominant role is played by the isoscalar 1"^ channel, which couples protons and neutrons in the same 
orbitals and in spin-orbit partners. As a consequence we find that, for ^'^'^^Fe, the proton pair matrix in the f^/2 
orbital, A/33, is larger than in the P3/2 orbital, although the latter is favored in energy. For ^^Fe, this ordering is 
inverted, as the increasing number of neutrons in the p^/2 orbital increases the proton pairing in that orbital. 

After diagonalization of M, the J = proton pairing strength is essentially found in one large eigenvalue. Fur- 
thermore, we observe that this eigenvalue is significantly larger than the largest eigenvalue on the mean-field level 
(Fermi gas), supporting the existence of a proton pair condensate in the ground state of these nuclei. The situation 
is somewhat different for neutrons. For ^^Fe, only little additional coherence is found beyond the mean-field value, 
reflecting the closed-subshell neutron structure. For the two other isotopes, the neutron pairing exhibits two large 
eigenvalues. Although the larger one exceeds the mean-field value and signals noticeable additional coherence across 
the subshells, the existence of a second coherent eigenvalue shows the shortcomings of the BCS-like pairing picture. 

It has long been anticipated that J = 0^ proton-neutron correlations play an important role in the ground states 
of = Z nuclei. To expl ore these correlations, w e have performed SMMC calculations of the N = Z nuclei in the 
mass region A = 48 — 56 (Langanke et ai, 1997a). Note that for these nuclei the pair matrix in all three isovector 



0+ channels essentially exhibits only one large eigenvalue, related to the fj/2 orbital. We will use this eigenvalue as 
a convenient measure of the pairing strength. As the even-even N = Z nuclei have isospin T = 0, the expectation 
values of A''' A are identical in all three isovector 0+ pairing channels. This symmetry does not hold for the odd-odd 
N — Z nuclei in this mass range, which have T = 1 ground states, and {A'' A) can be different for proton-neutron 
pairs than for like-nucleons pair (the expectation values for proton pairs and neutron pairs are identical). We find the 
proton-neutron pairing strength significantly larger for odd-odd N = Z nuclei than in even-even nuclei, while the 0"*" 
proton and neutron pairing shows the opposite behavior, in both cases leading to an odd-even staggering, as displayed 
m Fig. |l|. This stag gering is caused by a constructive interference of the isotensor and isoscalar parts of A'' A in the 
odd-odd N = Z nuclei, while they interfere destructively in the even-even nuclei. The isoscalar part is related to the 
pairing energy, and is found to be about con stant for the nuclei studied here. Similar behavior was also demonstrated 

This model is analytic, but shows the same 



19961 , |1998| ) 



in a simplified 5*0(5) seniority-like model (Engcl et al. 
trends as the shell model results. Due to other correlations present in the shell model, such as the inclusion of several 
orbits, isoscalar pairing, spin-orbit splitting, long-range correlations, deformation, etc., the shell model results are 
reduced in comparison to the simplified model. Pair i ng co r relations have also be en studied in heavier systems that 
require the presence of the Ogg/2 orbital (Dean et ai, 1997; Petrovici et at, 2000). In heavier odd-odd N = Z nuclei 
the g round state becomes a T = 1 (rather than T = 0), as was found experimentally in ''^Rb (Rudolph and et a], 
1996 ). The lowest T = and T — 1 states in these systems are very close in energy. Recently, mean-fie l d calc ulations 
that include both T = and T = 1 pairing correlations in odd-odd N = Z nuclei ( ^atula and Wysi , 2001) showed 
that the interplay between quasiparticle excitations (relevant for the case of T = states) and isorotations (relevant 
for the case of T = 1 states) explains the near degeneracy of these states. 
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D. Proton-neutron pairing and the Wigner energy 



So far, the strongest evidence for np pairing comes from the masses of N=Z nuclei. An additional binding (the 
so-called Wigner energy) found in these nuclei manifests itself as a spike in the isobaric mass parabola as a function 
of r^=i(A^ - Z) (see the review ( ^elde4 |1996D and references quoted therein). Gross estimates of the magnitude 
of the Wigner energy come from a large-scale fit to experimental binding energies with the macroscopic-microscopic 
appr oach (Krappe et al., 1979; Myers and Swiatecki, 1966) and from the analysis of experimental masses (Jensen et al. 



1984). Several early attempts were made to incorporate the effect of neutron-proton pair correlations in light nuclei 
in qu asiparticle theory (for an early review see ( poodmaij , |1979| )) with varying success. Satula et al, (patula et al. 



1997) presented a technique to extract the Wigner energy directly from the experimental data and gave empirical 
arguments that this energy originates primarily from the T=0 part of the effective interaction. To obtain deeper 
insight into the structure of the Wigner term, they applied the nuclear shell model to nuclei from the sd and fp shells. 

The Wigner energy, E^y, is believed to represent the energy of collective np-pairing correlations. It enters the 
semi-empirical mass formula (see e.g. (Krappe et al, 1979| )) as an additional binding due to the np-pair correlations. 
The Wigner energy can be decomposed into two parts: 



Ew = W{A)\N -Z\+ d{A)TTnp5 



npONZ, 



(64) 



The — Z|-dependence in Eq. ( |64[ ) was first introduced by Wigner ( Wigner , 1937 ) in his analysis of the SU(4) 
spin-isospin symmetry of nuclear forces. In the supermultiplet approximation, there appears a term in the nuclear 
mass formula which is proportional to Tgs(Tgs -I- 4), where Tgs denotes the isospin of the ground state. Empirically, 



Tgs—\Tz\ for most nuclei except for heavy odd-odd N=Z systems (Janecke, 1965| ; Zeldes and Liran, 1976). Although 
the experimental data in dicate th a t the SU(4) symmetry is severely broken, and the masses behave according to the 
Tgs (Tgs -1-1) dependence ( Janecke , 1965 ; Jensen et al. , 1984 ), the expression of Eq. ( |64|) for the Wigner energy is still 
very useful. In particular, it accounts for a non-analytic behavior of nuclear masses when an isobaric chain crosses the 
N=Z line. An additional contribution to the Wigner term, the d-term in Eq. (|6^), represents a correction for N=Z 
odd-odd nucl ei. Th e oretical justi f icatio n of Eq. (M) has been given in terms of basic properties of effective shell-model 
i nteractions ( Talmi, 1962; Zelde^ |1996 ), and also by using simple arguments based on the number of valence np-pairs 
(lensen et al, 1984; Myer4 1977). The estimates based on th e number of np pairs in identical spatial orbits suggest 



1984): d/W^=0.56±0.27. 



that the ratio d/W is constant and equal to one (Myers, 1977). A different estimate has been given in ( Jensen et al. 



In the work of Satula ( Satula et al. , 1997), the Wigner energy coefficient W in an even-even nucleus Z=N—^ was 
extracted by means of the indicator: 



W{A) = 6Vn: 
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The d-term in an odd-odd nucleus, Z^N^q [Eq. (g)], can be extracted using another indicator: 
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where the double-difference formula from (Zhang et al., 198!;) is 



5Vnp{N,Z) = ^{B{N,Z)-B{N-2,Z)-B{N,Z-2) + B{N-2,Z-2)} 

_ d^B 
dNdZ' 



(67) 



Although the recipe for these third-order mass difference indicators is not unique, the results appear to be very weakly 
dependent on the particular prescription used. 

To visualize the influence of the T=0 part of the effective nuclear interaction on the Wigner term, we have performed 
a set of shell-model calculations while switching oj9^ sequentially the J = 1, 2, Jmax, T=0 two-body matrix elements 
Uij2JT\H\j[j2JT) of the shell-model Hamiltonian H for different values of Jmax- Figure ^shows a ratio ew/^^^^^, 
where e^*^' denotes the result of full shell- model calculations versus Jmax- The calculations were performed for two 
representative examples, namely, the /p-shell nucleus *^Cr and the sd-shell nucleus ^"^Mg. The largest contribution to 
the Wigner energy comes from the part of the T=0 interaction between deuteron-like (J=l) and 'stretched' [J— 5 {sd) 
and 7 (pf)] pairs. The importance of these matrix elements is well known; it is prec isely for J~l and stretched pair- 
states that experimentally determined effective np T=0 interactions are strongest (Anataraman and Schiffer, 1971 
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Molinari et ai, 1975; Schiffei, 1971). Note also that the deuteron-Hke correlations contribute more strongly to ew in 



sd-nuclei than in /p-nuclei, and that matrix elements corresponding to intermediate values of J give non-negligible 
contributions. This reveals the complex structure of the Wigner energy and suggests that models which ignore high- J 



components of the np interaction (e.g., by considering only J=0, T=l and J=l, T~0 np pairs ( Evans et al. 
are not too useful for discussing the actual np pair correlations. 



1981D ) 



E. Thermal properties of p/-shell nuclei 



The proper ti es of nuclei at finite temperatur es are of c onsid e rable experiment a l (for reviews, see (jSnover, 1986 



Suraud et ai , 1989| )) and theoretical interest ( Alhassid , 1991; Egido and Ring , 1993 ). How thermal excitations 



influence the pairing properties will be the main focus of this section. We address this topic in Sec. ^ as well but 
with an emphasis on an analysis based on experimental data on level densities. 

Theoretical studies of nuclei at finite temperature have been based mainly on mean-field approaches and thus only 
consider the temperature dependence of the most probable configuration in a given system. These approaches have 
been criticized due to their neglect of quantum and statistical fluctuations (Dukclsky et ai, 1991). The SMMC 



method does not suffer this defect and allows the consideration of model spaces large enough to account for the 
relevant nucleon-nucleon correlations at low and moderate temperatures. 

SMMC calculations were pe rformed to s tudy the thermal propcirties of several even-even and odd-A nuclei in 



the mass region A = 50 — 60 ( Dean et al , 1994; Langanke et al , 1996 ) in an /p-shell model space using realistic 



interactions. More recently, Alhassid et al, carried out thermal calculations in a larger model-space which included 
the OgQ/2 shell. As a typical example, we discuss in the following our SMMC results for the nucleus ^'*Fe, which is 
very abundant in the presupernova core of a massive star. 

Our calculations include the complete set of 1753/2,1/20/7/2,5/2 states interacting through the realistic Brown-Richter 



1991) 



Hamil tonian (Richter et al 
I98I4 H) give essentially the same results 



(SMMC calculations using the modified KB3 interaction (Poves and Zukci 
Some 5 x 10^ configurations of the 8 valence neutrons and 6 valence 
protons moving in these 20 orbitals are involved in the canonical ensemble. The results presented below have been 
obtained with a time step of A/3 = 1/32 MeV~^ using 5000-9000 independent Monte Carlo samples. 

The calculated temperature dependence of various observables is shown in Fig. |2^. In accord with gener al thermo- 
dynamic principles, the internal energy U steadily increases with increasing temperature (Dean et al, 1994). It shows 
an inflection point around T w 1.1 MeV, leading to a peak in the heat capacity, C = dU / dT , whose physical origin we 
will discuss below. The decrease in C for T > 1.4 MeV is due to our finite model space (the Schottky effect ( pivitare"se 
et a[], [1989 )); we estimate that limitation of the model space to only the p/-shell renders our calculations of ^**Fe 
quantitatively unreliable for temperatures above this value (internal energies U ^ lb MeV). The same behavior is 
apparent in the level density parameter, a = C/2T. The empirical value for a is A/^ MeV — 6.8 MeV~^, which is in 
good agreement with our results for T « 1.1-1.5 MeV. 



More recent calculations confirm these basic findings. Liu and Alhassid calculated ( Liu and Alhassid , 2001) the 
heat capacity for iron isotopes in a complete 0/lp — gg/2 model space. They used a phenomenological pairing-plus- 
quadrupole model for the two-body interaction and found that the pairing transition in the heat capacity is correlated 
with the suppression of the number of spin-zero neutron pairs as the temperature increases. The results were obtained 
using a novel method to calculate the heat capacity that decrease the statistical error bars in t he calculation. We 
show results of this calculation for Fe isotopes in Fig. |2l|. While the original calculations of ( Dean et al\ 1994) 
indicate a possible phase transition (along with a pairing collapse in the measure d (A^A) pairing ex pectation), this 
effect appears to be delayed to more neutron- rich nuclei in the calculations of (Liu and Alhassid, ^001 ). Several 
factors likely contribute to this difference. First, the interactions are obviously different. Second, the extrapolation 
techniques used for realistic interactions likely over-estimate the influence of pairing in the region between 0.5 and 
1.0 MeV. Finally, the model space is smaller in the early calculation, although the Schottky pea k is seen to appea r 
at about 1.4 MeV. This makes the interpretation of the low temperature peak more difficult in ( Dean et al , 1994). 
Neverthe less, it should be clear that both the original calculations with realistic interactions and the more recent 
work in ( Liu and Alhassid , 2001) both indicate interesting physics related to pairing phenomena in the T = 1.0 MeV 
region. 

We also show in Fig. ^ the expectation values of the BCS-like proton-proton and neutron-neutron pairing fields, 
(At A). At low temperatures, the pairing fields are significantly larger than those calculated for a non- interacting 
Fermi gas, indicating a strong coherence in the ground state. With increasing temperature, the pairing fields decrease, 
and both approach the Fermi gas values for T « 1.5 MeV and follow it closely for even higher temperatures. Associated 
with the breaking of pairs is a dramatic increase in the moment of inertia, / = (J^)/3T, for T = 1.0-1.5 MeV; this 
is analogous to the rapid increase in magnetic susceptibility in a superconductor. At temperatures above 1.5 MeV, / 
is in agreement with the rigid rotor value, 10.7?i^/MeV; at even higher temperatures, it decreases linearly due to our 
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finite model space. 

Although the results discussed above are typical for even-even nuclei in this mass region (including the N = Z 
nucleus ^^Fe), they are not for odd-odd N = Z nuclei. This is illustrated in Fig. ^ which shows the thermal behavior 
of several observa bles for ^"Mn (N ^ Z = 25), calculated in a SMMC study within the complete p/-shell using the 
KB3 interaction (Poves and Zukci, 198l4 |b|)- A closer inspection of the isovector J — and isoscalar J = 1 pairing 



correlations holds the key to the understanding of these differences. The J = isovector correlations are studied 
using the BCS pair operators, Eq. ( |5^ ) with a similar definition for proton-neutron pairing. For the isoscalar J = 1 
correlations, we have interpreted the trace of the pair matrix M'^^^ (defined in Eq. (^0|)) as an overall measure for 
the pairing strength, 

PL = Y.>'i^Y.Mic.- (68) 

/3 a 

Note that at the level of the non-interacting Fermi gas, proton-proton, neutron-neutron, and proton-neutron J — Q 
correlations are identical for N = Z nuclei. However, the residual interaction breaks the symmetry between like-pair 
correlations and proton-neutron correlations in odd-odd N = Z nuclei. As is obvious from Fig. at low temperatures 
proton-neutron pairing dominates in ^"^Mn, while pairing among like nucleons shows only a small excess over the Fermi 
gas values, in strong contrast to even-even nuclei. 

A striking feature of Fig. ^ is that the isovector proton-neutron correlations decrease strongly with temperature 
and have essentially vanished at T = 1 MeV, while the isoscalar pairing strength remains about constant in this 
temperature region (as it does in even-even nuclei) and greatly exceeds the Fermi gas values. We also note that 
the pairing between like nucleons is roughly constant at T < 1 MeV. The change of importance between isovector 
and isoscalar proton-neutron correlations with temperature is nicely reflected in the isospin expectation value, which 
decreases from < >= 2 at temperatures around 0.5 MeV, corresponding to the dominance of isovector correlations, 
to < >= at temperature T = 1 MeV, when isoscalar proton- neutron correlations are most important. 

The temperature dependence of the excitation energy E = (H) in the odd-odd nucleus ^°Mn is significantly different 
than that in even-even nuclei. The difference is due to the uniqueness of isospin properties in odd-odd N ~ Z nuclei. 
It is only in these nuclei that one finds states of different isospin, T = 1 and T = 0, that are close to each other at low 
excitation energies. The ^°Mn ground state is T = 1, = 0, an d = 0+. In that state pn pair correlations dominate, 
and the like-particle correlations are reduced (Langanke et al, 1997a| ). However, at relatively low excitation energy 



these nuclei exhibit a multiplet of T = states with nonvanishing angular momenta. These states contribute efficiently 
to corresponding thermal averages. On the other hand, it follows from isospin symmetry that in the T = states all 
three pairing strengths (in Tz) must be equal. Thus, at temperatures where the T = states dominate the thermal 
average, the pn pair correlations are substantially reduced when compared to ground state values. This argument 



appears to be a generic feature of odd-odd N = Z nuclei beyond °Ca. For a further discussion see (Langanke et al. 
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1997b). For a different point of view from the perspective of mean-field calculations, see, for example (Ropke et al. 



F. Pair correlations and thermal response 

All SMMC calculations of even-even nuclei in the mass region A = 50 — 60 show that the BCS-like pairs break 
at temperatures around 1 MeV. Three observables exhibit a particularly interesting behavior at this phase transi- 
tion: a) the moment of inertia rises sharply; b) the Ml strength shows a sharp rise, but unquenches only partially; 
and c) the Gamow- Teller strength remains roughly constant (and strongly quenched). Note that the B{M1) and 
B[GT-\-) strengths unquench at temperatures larger than « 2.6 MeV and in the high-temperature limit approach the 
appropriate values for our adopted model space. 

( p^^anganke et al. , 199(^ ) has studied the pair correlations in the four nuclei ^^"^^Fe and ^^Cr for the various isovector 



and isoscalar pairs up to J = 4, tentatively interpreting the sum of the eigenvalues of the matrix A/'' ^ as an overall 
measure for the pairing strength. Note that the pairing strength, as defined in is non-zero at the mean-field level. 
The physically relevant pair correlations Pc'orr ^^'S then defined as the difference of the SMMC and mean-field pairing 
strengths. 

Detailed calculations of the pair correlations have been performed for selected temperatures in the region between 
T — 0.5 MeV and 8 MeV. Fig. p3| shows the temperature dependence of those pair correlations that play an important 
role in understanding the thermal behavior of the moment of inertia and the total Ml and Gamow- Teller strengths. 

The most interesting behavior is found in the J — Q proton and neutron pairs. There is a large excess of this 
pairing at low temperatures, refiecting the ground state coherence of even-even nuclei. With increasing temperature, 
this excess diminishes and vanishes at around T = 1.2 MeV. We observe further from Fig. ^ that the temperature 
dependence of the J = proton-pair correlations are remarkably independent of the nucleus, while the neutron pair 
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correlations show interesting differences. At first, tfie correlation excess is smaller in the semimagic nucleus ^^Fe than 
in the others. When comparing the iron isotopes, the vanishing of the neutron J — correlations occurs at higher 
temperatures with increasing neutron number. 

The vanishing of the J — proton and neutron pair correlations is accompanied by an increase in the correlations 
of the other pairs. For example, the isovector J — proton-neutron correlations increase by about a factor of 3 after 
the J = proton and neutron pairs have vanished. The correlation peak is reached at higher temperatures with 
increasing neutron number, while the peak height decreases with neutron excess. 

The isoscalar proton-neutron J = 1 pairs show an interesting temperature dependence. At low temperatures, when 
the nucleus is still dominated by the J = proton and neutron pairs, the isoscalar proton-neutron correlations show 
a noticeable excess but, more interestingly, they are roughly constant and do not directly reflect the vanishing of the 
J — proton and neutron pairs. However, at T > 1 MeV, where the J = proton- and neutron-pairs have broken, 
the isoscalar J — 1 pair correlations significantly increase and have their maximum at around 2 MeV, with peak 
values of about twice the correlation excess in the ground state. In contrast to the isovector J = proton-neutron 
pairs, the correlation peaks occur at lower temperatures with increasing neutron excess. We also observe that these 
correlations fade rather slowly with increasing temperature. 

A further discussion on thermal properties through recent experimental information on level densities will be given 
in Sec. 



IV. RANDOM INTERACTIONS AND PAIRING 



We have seen that all even-even nuclei have a J'^ = 0^ ground state. Pairing in even-even systems was also 
shown in previous sections to be a major contributor to the ground-state correlations. Furthermore, a property like 
the O"*" nature of all even nuclei can be explained within the simple seniority model, based itself on the short-range 
nature of the effective interaction. It is therefore interesting to see whether such a general property is specific of this 
hamiltonian or whether it could also emerge from a ra ndom ensemble of ro tationally and isospin invariant random 
two-body interactions. This question was first posed in ( lohnson et ai , 1998 ), where the low- lying spectral properties 
of random interactions were first studied from the shell-model perspective. Several interesting results were obtained 
including highly likely 0"*" ground state s emerging from the random ensembles, a n enhanced phonon coll ectivity, 
strongly correlated pairing phenomena ( [lohnson et ai , 2000), odd-even staggering (Papenbrock et al, 2002), and a 
likelihood of generating rotational and vibrational spectra. 

Similar results were also obtained in the interacting boson model (IBM) (Bijker and Frank, 200C). The IBM 
Hamiltonian used in these calculations is given by 



H 



ehd - nQix) X Qix) 



(2) 

,(X) = (stj+dt^) +x[d^d^ 



(69) 
(70) 



where only spin (s, monopole) and spin 2 {d, quadrupole) bosons are considered. This interaction is randomized 
by introducing a scaling parameter rj — e/ [e + 4k(7V — 1)] and x = 2x/'\/7, and choosing x and rj randomly on the 
intervals [—1,1] and [0, 1], respectively. The IBM calculations also gave a predominance of 0+ ground states as well as 
strong evidence for the occurence of both vibrational and rotational band structures. Within the shell model, these 
structures appear within a continuum of rotational bands, but the nature of the IBM model restricts the structures 
to be of these two particular forms. In this section we will briefly review the present status of research into this 
interesting phenomena. 

For fermions, we define the two-body matrix elements VjT{o.b, cd) through an ensemble of two-particle Hamiltonians 
and require that the ensemble be invariant under changes in the basis of two-particle states. This is achieved by taking 
the matrix elements to be Gaussian distributed about zero with the widths possibly depending on J and T such that: 



{Via') ^CJ^,TAI + Sa^a')v\ 



(71) 



Here v is an overall energy scale that we generally ignore (except for scaling single-particle energies for the RQE-SPE 
defined below). The coefficients cj then define the ensemble. We emphasize that J, T refer to quantum numbers of 
two-body states and not of the final many-body states (typically 4-10 valence particles). 

Several basic ensembles may be defined by the choice of the form of the cj^t coefficients and the single- 
particle Hamiltonian, if present. One ensemble is called the Random Quasiparticle Ensemble (RQE). In this case 
cj,T = [{'2T + 1)(2J-|- 1)]~^. This relation between the cj.t, which was discussed in ( polmson et al\ , 1998), came 
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from imposing on the ensemble the constraint that it should remain the same for the particle-particle interaction as 
for the particle-hole interaction. A different ensemble in this class is called the two-body random ensemble (TBRE) 
for which cjt = constant. Historically, this wa s the first two-part i cle ra ndom ensemble to be employed in studying 
statistical properties of many-particle spectra ( French and Wong , 1970 ). These two ensembles assume degenerate 
single-particle energies. Realistic interactions do have nondegenerate single-particle energies that will, in principle, 
affect variou s spectral p r opertie s. For calculations in the sd shell one uses single-particle energies from the Wildenthal 
interaction (Wildenthal, 1984a), scaling v = 3.84 MeV to best match the widths of the two-particle matrix elements. 
The resulting interactions with the single-particle splitting included are called the RQE-SPE and TBRE-SPE. 

The first, and perhaps most striking, feature of all of these random interactions is the preponderance of = 0-|- 
ground states. In Fig. |2j we show the distribution of ground-state spins in the various ensembles for the two systems 
^•^Ne and ^^Mg. We generated 1000 random interactions from each ensemble. These results are typical and consistent 
with calculations with only one type of particle (for example, neutrons only), or fermions in which the I ■ s force is 
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not present (Kaplan et al. 

Fig. |2j also shows that the even spins are preferred. In some cases, higher even spin states are preferred over 
medium spin states. For example, in ^"^Mg, the S 
to lower slightly the number of 0+ ground states 
in each case. 

Various research efforts have been undertaken to understand the preponderance of the 0^ ground-state. Ensembles 
of interactions derived from a Gaussian unitary ensemble (GUE) distribution are not time-re versal invariant, bu t both 
the Gaussian orthogonal ensemble (GOE ) and GUE random ensembl es yield 0+ dominance (Bijker et ai, 1999). This 



state is preferred over the 6"*". The single-particle splitting tends 
The RQE clearly obtains the highest number of 0+ ground states 



apparent paradox was recently resolved ( Velazqucs and Zuker , 2002| ) by noting tha t the operator commutes with 
the T time-reversal operator for either ensemble. For bosons, Kusnezov ( [Kusnezov , 200C) was able to map the C/(4) 
vibron model onto random polynomials on the unit interval. Kusnezov was then able to show analytically the origin 
of 0+ ground states. While the C/(4) model is extremely simplified and only describes bosons (rather than fermions), 
it points to the interesting link between random polynomials and the two-body interaction. The 0+ predominance in 
the fermion case was recently studied by Mulhall et al (Mulhall et ai, 2000). These authors used a single j-shell to 



show that statistical correlations of fermions in a finite system with random interactions drive the ground state spin 
to its minimum or maximum value. The effect is universal and related to the geometric chaoticity, or the assumption 



of pseudorandom coupling of angular momentum (Zclevinsky et al, 199(:), of the spin coupling of individual fermions 



While a rigorous derivation of these findings for general orbital schemes is not yet available, the research is pointing 
towards an understanding of why an ensemble of random interactions posesses predominantly 0^ ground states. 

A second feature concerns the likelihood of finding rotational or vibrational spectra from the ensembles. The 
relevant measure for these states is the ratio of the first 4+ excitation energy to the first 2+ energy. This ratio, 
p — [£'(4+) — E{0~^)]/[E{2^) — £'(0+)] is 2 for a vibrational spectrum and 10/3 for a rotational one. The results from 
the RQE show a broad distribution of various kinds of spectra peaked towards vibrational spectra. Using random 
interactions in the IBM, virtually all random interactions yield a vibrational or rotational spectrum in nearly equal 
proportions. The difference is due to the restricted nature of the random I BM interaction in whi ch only s and d 
boson pairs and couplings are included in the Hamiltonian. Kusnezov et al. (Kusnezov et al, 200C) confronted the 



results of the IBM model with known experimental data and found two interesting results. They found that both the 
interaction and the number of relevant valence nucleons were key to understanding the distribution of pibm- They 
also found that experimental data favor rotational spectra over vibrational spectra and were able to place limits on 
the choice of random variables that would allow for a reproduction of experimental data. 

Of the three general properties of random interactions we dis cuss here, the enh ancement of the B{E2) strength is not 
spontaneously produced by our choices of rand om interactions ( Horoi et al. , 2001). This is particularly true for strongly 
deformed nuclei. The problem can be cured ( Velazques and Zukei , 2002) by choosing a constant displacement of all 
the matrix elements which is essentially the same as adding some coherence to the choice of the random Hamiltonian. 
Velazques and Zuker ( Velazques and Zuker , ^002 ) were able to do this and showed how one may obtain good rotational 
spectra from the displaced-TBRE. 

A third feature of the random shell-model interactions involves the pair content of the 0^ ground-state wave func- 
tions. The pairing content of the wave function was measured by calculating the pair-transfer operator. Interestingly, 
for a given interaction, the same coherent pair connected several even-even nuclei in a given isotopic chain. This 
feature appears to be robust. On the other hand, studies (Mulhall et al., 2000) made in a single j shell relate the 



origin of regularities in the spectra to incoherent interactions rather than to coherent pairing. The origin of order in 
the spectra, attributed to geometric chaocity, tends to prove that the role of pairing is minimal. In order to better 
understand these two seemingly confli cting observations, a further analys is of the pairing properties of the system has 
been performed by Bennaceur et al. (Bennaceur et al., submitted, 2002) who compared shell-model results to those 
obtained from Hartree-Fock-Bogoliubov (HFB) calculations using the same set of random interactions. The aim of 
this study was to determine whether the interactions support static pairing or whether the effect is more dynamical. 
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Because the HFB solution generally breaks all the symmetries required by the many-body Hamiltonian, it is not a 
physical state, but an indicator of the intrinsic structure in the many-body system. 

In the HFB approach one does not explore the full Hilbert spa ce; the trial funct i on is constrained to be a superpo- 



sition of pairs of single particles. Moreover, the Hamiltonian of ( Bennaceur et al , submitted, 2002 ) omits the terms 
which represent the residual interaction between quasiparticles, and proton-neutron pairing is not taken into account. 
Finally, the HFB approach is an unprojected variational method, so J is not a good quantum number and neither are 
the particle numbers N and Z . In the shell model the particle numbers are well defined, while in HFB approximation, 
only the average particle numbers are constrained. 

The HFB approximation describes only static pairing. On the contrary, in the SM model picture, the wave functions 
contain all the possible correlations inside a given model space. For that reason, a shell-model ground-state wave 
function can show some strong pairing properties while the corresponding HFB solution can be totally unpaired. In 
that case the pair structure can be due to dynamical pairing which cannot be described in the HFB method but 
require going beyond the mean field. 

The three systems ^''Mg, ^^Ne, and ^''O were considered. Only those interactions generated from the RQE-SPE 
ensemble that lead to a SM ground state with J'^ = 0+ were used, and the pairing properties of the ground state 
wave functions were investigated. The simultaneous use of the shell model and the HFB approximation gives a better 
understanding of the pairing induced by the random interactions. 

In order to investigate the pairing properties of the shell- model solutions, one introduces the pair transfer coefficient 

< P+ > ^ < A\P+\A -2>=Y^<A\ [alat]l\A - 2 > , (72) 

where \A > and \A — 2> represent the ground states of the A and A— 2 particle systems obtained from the shell model 
(isospin T = 1 is understood). This quantity is compared to the mean pairing strength in the HFB approximation 
< K > defined by 

<K>=TrK. (73) 

The three systems considered correspond in the (sd) shell to N — Z ~ 0, N — Z = 2, and N — Z — 4. 

If the dominance of J'" = 0+ ground state is due to the pairing properties of the system, then we expect to obtain 
a significant value for < k > in most of the cases. Moreover if the pairing plays an important role for the structure 
of the ground state, then it must be related with the pair transfer coefficient, and in that case one can also expect a 
clear correlation between < k > and < P+ > . 

In Fig. E3, we represent the distribution of the number of interactions according to the results obtained for < P+ > , 



< K>, and Q2 (Bennaceur et al, submitted, 2002). In the case with only one kind of active particle in the model 
space {N — Z = 4), almost all the interactions lead to a strong pair transfer coefficient. This property can be related 
to the fact that when we consider only one kind of particle, the deformation effects are (almost) zero, and pairing 
can develop more easily. Netherless, for the three sets, we see once again that <P+ > has a significant value in most 
cases. The distribution of the number of interactions according to the static pairing strength, measured via < k>, 
follow the same evolution (i.e. the number of interactions that give a significant value of < k > increases with N ^ Z), 
but the number of interactions for which < k > is small is always important. It is then very unlikely to relate the 
origin of the spin of the ground states with the static pairing created by the interactions. 

It is also instructive to consider the evolution of the plots when one changes the asymmetry of the system from 
N — Z = AtoN— Z = t). In the mean-field description (center and lower part of Fig. ^), the pairing strength is 
concentrated into regions < k > ~ and <k>~ 0.5 to 2. This property does not strongly change as a function of 
N ~ Z. Nevertheless, we notice that < k > is more often close to zero and the non-zero values are less scattered when 
N = Z . This effect can probably be attributed to the deformations which play a more important role when N ^ Z . 
When the pairing is weak, the deformation effects prevail and so decrease the pairing strength in the region between 
and 1.5. 

In the shell-model description (upper part of Fig. ^5|), we notice a clear evolution of the pair transfer coefficient 
with the asymmetry of the system. This property seems to be mainly due to the deformation effects. Indeed for the 
system with TV — Z = 4, the coefficient < P+ > is peaked at around 2.5, and no interactions give a value close to zero. 
The cases with N — Z = 2 and N — Z — are similar and indicate that the T = part of the pairing interaction 
does not play a crucial role in these systems. This last remark tends to show that the T — part of the pairing 
interaction, like the T = 1 part, does not play a significant role for the r elative abundance of 0+ ground states. This 



conclusion is in agreement with the statement made by Mulhall et al. (Mulhall et al, 2000), in which the origin of 
the abundance of the 0"*" states in the even fermion systems is related with the geometric chaocity rather than with 
the pairing properties of the system. 
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V. THERMODYNAMIC PROPERTIES AND PAIRING CORRELATIONS IN NUCLEI 
A. Introduction 



One of the most interesting problems in the context of phase transitions in small systems is the question of the 
existence and classification of a possible phase transition from a hadronic phase to a quark-gluon plasma in high- 
energy physics. The answer to this question has far-reaching consequences into many other fields of research like, 
e.g., cosmology, since it has been argued that hadronization of the quark-gluon plasma should be a first-order phase 
transition in order to allow for possible super cooling and consequently the emergence of large-scale inhomogeneities 
in the cosmos within the inflationary big-bang model. 

In nuclear physics, different phase transitions have been discussed in the literature. A first-order phase transition has 
been reported in the multifragmentation of nuclei (D'Agostino et ai, 2000), thought to be the analogous phenomenon 



in a finite system to a liquid-gas phase transition in the thcrmodynamical limit. A pivotal role in these studies is 
playe d by the presence of a convex intruder in the microcanonical entropy curve (Gross, 1997; [Gross and Votyakov , 
2000| ). This leads to a negative branch of the microcanonical heat capacity which is used as an indicator of a first-order 
phase transition in small systems. Negative heat capacities have indeed been observed in the multifragmentation of 
atomic nuclei, thoug h the heat capacity curve has not been derive d directly from the caloric curve, but by means of 
energy fluctuations (Chomaz et al., 2000; D'Agostino et al\ , 20001). Another fi nding of a negative branch of the heat 
capacity curve has been in sodium clusters of 147 atoms ( [M. Schmidt et ai , 2001), indicating a possible first-order 
phase transition. On the other hand, it is not clear whether the observed negative heat capacit ies are simply not du e 



to the changing volume of the system under study that is progressively evaporating particles (Moretto et ai, 2001) 



In general, great care should be taken in the proper extraction of temperatures and other thcrmodynamical quantities 
of a multifragmenting system. 

Another transition discussed for atomic nuclei has been anticipated for the transit ion from a phase with strong 
pairing correlations to a phase with weak pairing correlations ( ^ano and Yamasaki , 1963). Early schematic calculations 
have shown that pairing correlations can be quenched by temperature as well as by the Coriolis force in rapidly 
rotating nuclei ( poodman , 1981a b; Shimizu et ai, 1989| ; [Tanabe and Sugawara- Tanabt^ , [1980| , |1982| ). This makes the 
quenching of pairing correlations in atomic nuclei very similar to the breakdown of superfiuidity in ^He (due to rapid 
rotation and/or temperature) or of superconductivity (due to external magnetic fields and/or temperature). Recently, 
structures in the heat capacity curve related to the quenching of pairing correlations have been obtained within 
t he relativis t ic me an-field theory (Agrawal et al, 200C| , |2001| ), the finite-ter nperature randoni phase approximation 
(Dinh Dang 



190C), the finite-tem perature Hartree-Fock-Bogoliubov theory ( Egido et al , 2000|), and the s hell- model 



Mont e Car l o (SMMC) approa ch ( [Dean et al] , |1995| ; |Liu and Alhassi4 |2001| ; [Nakada and Alhassi4 |1997| ; [Rombout, 
et a|] , 199S ; [White et al , 2000 ) . An S-shaped structure in the heat capacity curve could also be observed experimentally 
( Schiller et al , 2001) and has been interpreted as a fingerprint of a second-order phase transition in the thermodynamic 
limit from a phase with strong pairing correlations to one with no pairing correlations. For finite systems there will be 
a gradual transition from strong pairing correlations to weak pairing correlations, implying a finite order parameter, 
here the pairing gap A. Indeed, the analogy of the quenching of pairing correlations in atomic nuclei with the 
breakdown of superfluidity in '^He and the breakdown of superconductivity suggests a second-order phase transition 
and a schematic calculation might support this assumption (see the discussion below). Interestingly, similar structures 
of the heat capacity curve as observed for atomic nuclei in ( ^chiller et al\ , |200l| ) have been seen in small metallic grains 



undergoing a second-order phase transition from a superconductive to a n ormal conductive phase (Black et al, 1996 
1997; von Delft and Ralph, ^001 ; Lauritzen et al, 1993; Ralph et al, 1995[ ), thereby supporting the analogous findings 
for atomic nuclei. On the other hand, breaking of nucleon pairs has been experimentally shown to cause a series of 
convex intruders in the microcanonical entropy curve of rare earth nuclei (Melby et al, 1999, 2001), leading to several 
negative branches of the microcanonical heat capacity. This finding might, in analogy to the discussion of nuclear 
multifragmentation, be taken as an indicator of several first-order transitions. Other possible phase transition-like 
behaviors are, e.g., shape transitions from a collective to an oblate aligned-particle structure at higher temperatures 
(see for example the recent work of Ma et al (Ma. et al, 2000| )). 



For a finite isolated many-body system such as a nucleus, the correct thcrmodynamical ensemble is the micro- 
canonical one. In this ensemble, the nuclear level density, the density of eigenstates of a nucleus at a given excitation 
energy, is the important quantity that should be used to describe thermodynamic properties of nuclei, such as the 
nuclear entropy, specific heat , and temper ature. Bethe first described the level density using a non-interacting fcrmi 
gas model for the nucleons ( Bethc[, 1936). Modifications t o this picture, such as the back-shifted fermi gas which 
includes pair and shell effects (Gilbert and Cameron, 1965) not present in Bethe's original formulation, are in wide 
use. We note that several approaches based on mean-field theory have recently been developed to investigate nuclear 
level densities including recent wor k that incorporates BCS pairin g into the mean-field picture to derive level densities 
for nuclei across the periodic table (Dcmctriou and Goriely, 2001). Other mean-field applications based on the Gogny 
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effective nucleon-nu cleon interact i on (w hich includes pairing due to the finite range of the interaction) have also been 



developed recently ( Hilaire et al , 2001 ). The level density^ p defines the partition function for the microcanonical 
ensemble and the entropy through the well-known relation S{E) — kBln{p{E)). Here ks is Boltzmann's constant 
and E is the energy. In the microcanonical ensemble, we could then extract expectation values for thermodynamical 
quantities like temperature T, or the heat capacity C. The temperature in the microcanonical ensemble is defined as 



(74) 



It is a function of the excitation energy, which is the relevant variable of interest in the microcanonical ensemble. 
However, since the extracted level density is given only at discrete energies, the calculation of expectation values like 
T, involving derivatives of the partition function, is not reliable unless a strong smoothing over energies is performed. 



This case is discussed in detail in (Mclby et ai, 1999) and below. Another possibility is to perform a transformation 
to the canonical ensemble. The partition function for the canonical ensemble is related to that of the microcanonical 
ensemble through a Laplace transform 



Z(/3) 



/•OO 

/ dEp{E)exp{-PE). 
Jo 



(75) 



Here we have defined /3 = l/fc^T. Since we will deal with discrete energies, the Laplace transform of Eq. (75) takes 
the form 



Z(/3)=^Ai?p(ii;)exp(-/3i?), 



(76) 



where AE is the energy bin used. With Z we can evaluate the entropy in the canonical ensemble using the definition 
of the free energy 



F{T) = -kBTlnZiT) = (E{T)) - TS{T). 



(77) 



Note that the temperature T is now the variable of interest and the energy E is given by the expectation value (E) as 
a function of T. Similarly, the entropy S is also a function of T. For finite systems, fluctuations in various expectation 
values can be large. In nuclear and solid state physics, thermal properties have mainly been studied in the canonical 
and grand-canonical ensemble. In order to obtain the level density, the inverse transformation 



-1 t^ioo 

p{E) = — dpZ{f3)exp{PE), 



(78) 



is normally used. Compared with Eq. (|75|), this transformation is rather difficult to perform since the integrand 
exp {f3E + In Z{[3)) is a rapidly varying function of the integration parameter. In order to obtain the density of states, 
approximations like the saddle-point method, viz., an expansion of the exponent in the integrand to second order 



1999; 


Koonin et al.. 


1997; 


White et al. 


, 2000) 



For the ideal Fermi gas (FG), this gives the following density of states 
exp {2\/aE) 



Pfg{E) 



SV48 



(79) 



A/B with dimension MeV ^, A being the mass number 



where a in nuclear physics is a factor typically of the order a 
of a given nucleus. 

To obtain an experimental level density is a rather hard task. Ideally, we would like an experiment to provide us 
with the level density as a function of excitation energy and thereby the 'full' partition function fo r the microcanonical 
ensemble. It is only rather recently that experimentalists have been able to develop methods (Henden et al., 1995; 
Tveter et al., [1996 ) for extracting level densities at low spin from measured 7-spectra. These measurements were 
perfor med at the O s lo Cy clotron Laboratory. Detection of gamma-rays are obtained with the CACTUS multidetector 
array ( Melby et al. , 1999 ) using the ('^He,a7) and ('^He,'^He') reactions on several rare-earth nuclei. Assuming that 



Hereafter we use p for the level density in the microcanonical ensemble. 
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the experimental analysis is correct, the resulting microcanonical level density reveals structures between 1-5 MeV of 
excitation energy and were interpreted as indications of pair breaking in these systems. 

The Oslo experimental results lead us to ask whether we can simultaneously understand the thermodynamic and 
pairing properties of a nuclear many-body system. We are also led to questions concerning the nature of phase tran- 
sitions in a finite many-body system. After these introductory words, we briefly delineate the aim of this section. In 
the next subsection we present experimental level densi ties f or several rare-earth nuclei together with a thcrmody- 
namical analysis and possible interpretations. In Subsec. V.C the simple pairing model of Eq. ( pT| ) is used in a similar 
analysis in order to see whether such a simplified pairing Hamiltonian can mimick some of the features seen from the 
experimental level densities. Since this is a simplified model, we also present results from shell-model Monte Carlo 
calculations of level densities in the rare-earth region with pairing-plus-quadrupole effective interactions in realistic 
model spaces. 



B. Level densities from experiment and thermal properties 



The Oslo cy clotron group has developed a me t hod t o extract nucl e ar level dens i ties at low spin fro m measured 



7-ray spectra (Hcndcn et al 





1995; Melby et al, 


1999 





Schiller et al, 2000a, 2001; Tveter et al, 1996). The main 



advantage of utilizing 7-rays as a probe for level density is that the nuclear system is likely thermalized prior to the 
7-emission. In addition, the method allows for the simultaneous extraction of level density and 7-strength function 
over a wide energy region. 

The experiments were carried out with 45 MeV ^He-projectiles accelerated by the MC-35 cyclotron at the University 
of Oslo. The experimental data were recorded with the CACTUS multidetector array using the (^He,a7) reaction on 
several rare-earth nuclei such as ^"^^Sm, ^^^Dy, ^^^Dy, -'^^^Er, -'^^^Yb, and ^^^Yb, yielding as result the nuclei ^'*^Sm, 
i6iDy, is^Dy^ leegj.^ ^ i^iyb, and ^^^yb. The ( ^He,^He') rea c tion wa s used to obtain the nuclei I'^^Sm and ^^^Er. For 
a critical discussion of the latter reaction see (3chiller et al, 2000b). The charged ejectiles were detected with eight 
particle telescopes placed at an angle of 45° relative to the beam direction. Each telescope comprises one Si Ai? front 
and one Si(Li) E back detector with thicknesses of 140 and 3000 /zm, respectively. 

From the reaction kinematics, the measured a-particle energy can be transformed to excitation energy E. Thus, 
each coincident 7-ray can be assigned to a 7-cascade originating from a specific excitation energy. The data are sorted 
into a matrix of {E, E^) energy pairs. The resulting matrix P{E, Ej), which describ es the primary 7-spectra obtai ned 
at initial excitation energy E, is factorized according to the Brink- Axel hypothesis (Axel, 1962; P.M. Brink, 1955) by 



P{E, E-y) cx p{E - E-y)a{E-y 



(80) 



where the level density p and the 7-energy-dependent function a arc unknown. The iterative procedure for obtaining 
t hese two functions and the assumptions behind the factorization of this expression are described in more detail in 
( Henden et al , 199?: ; Schiller et al , ^OOOa ). The experimental level density p{E) at excitation energy E is proportional 



to the number of levels accessible in 7-decay. For the present reactions, the spin distribution is centered around (J) 
4.4 h with a standard deviation of aj ^ 2.4 h. Hence, the entropy^ can be deduced within the microcanonical 
ensemble, using 



S{E) = kg \iiN{E) ^ kg In 



Po 



(81) 



where N is the number of levels in the energy bin at energy E. The normalization factor po can be determined from 
the ground state band in the even-even nuclei, where one has N{E) ^ 1 within a typical experimental energy bin of 
~ 0.1 MeV. 

The extracted entropies for the I6i,i62]-)y ^^^^ i71472y]-, nuclei are shown in Figs. p6| and p7| . In the transformation 
from level density to entropy, po was set to po ~ 3 MeV~^ in Eq. (^l]). The entropy curves are rather linear, but 
with small oscillations or bumps superimposed. The curves terminate around 1 MeV below their respective neutron 
binding energies due to the experimental cut excluding 7-rays with E^ < 1 MeV. All four curves reach S" ~ 13 fc^, 
which by extrapolation correspond to 5 '--^ 15 ks at the neutron binding energy i3„. 

Note that there is an entropy excess for the odd systems, since many low-lying states can be reached without 
needing to break a pair. The experimental level density can be used to determine the canonical partition function 



The experiment reveals the level density and not the state density. Thus, also the observed entropy reveals the number of levels. The 
state density can be estimated by pstate ~ (2J + l)pievel ~ 9-8 Plevel- 
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Z{T). However, in the evaluation of Eg. ([76|), one needs to extrapolate the experimental p curve to ^ 40 MeV. The 



back-shifte d level density formu la of ( von Egidy et al , 1988 ; Gilbert and Cameron , 1965 ) was employed (for further 
details see (Schiller et ai, 2000a)). From this semi-experimental partition function, the entropy can be determined from 
Eq. (|8^). The results are shown in Fig. |2^. The entropy curves show a splitting at temperatures below fc^T = 0.5 — 0.6 
MeV which reflects the experimental splitting shown in the microcanonical plots of Figs. ^ and p7| . 

The merging together of the entropy curves at roughly ksT = 0.5 — 0.6 MeV can also be seen in the analysis of 
the heat capacity in the canonical ensemble. The extract i on of the microcanonical heat capacity Cv{E) gives large 
fluctuations which are difficult to interpret (Melby et ai, 1999). Therefore, the heat capacity Cv{T) is calculated 



within the canonical ensemble as function of temperature T. The heat capacity is then given by 



Cv{T) = 



0{E) 
dT ' 



(82) 



The deduced heat capacities for the -'^^-'^'-'^^^Dy and ^'^^'^'^^Yh nuclei are shown in Fig. ^ All four nuclei exhibit similarly 
S-shaped Cy(T)-curves with a local maximum relative to the Fermi gas estimate at Tc ~ 0.5 MeV. The S-shaped curve 
is interpreted as a fingerprint of a phase transition-like behavior in a finite system from a phase with strong pairing 
correlations to a phase without such correlations. Due to the strong smoothing introduced by the transformation to 
the canonical ensemble, we do not expect to see discrete transitions between the various quasiparticle regimes, but 
only the transition where all pairing correlations are quenched as a whole. In the right panels of Fig. we see that 
Cv{{E)) has an excess in the heat capacity distributed over a broad region of ex c itatio n energy and is not giving a 



19991 ). 



clear signal for quenching of pairing correlations at a certain energy (Melby et 

In passing, we note that the results displayed in Fig. |2^ are similar to those of Liu and Alhassid ( Liu and Alhassid , 
200 1|) shown in Fig. |2~ 



C. Thermodynamics of a simple pairing model 

In this section, we will try to analyze the results from the previous subsection in terms of the simple pairing model 
presented in Eq. (|2l]). As stated in Sec. ||, seniority is a good quantum number, which means that we can subdivide 
the full eigenvalue problem into minor blocks with given seniority and diagonalize these separately. If we consider an 
even system of TV = 12 particles which are distributed over L = 12 two-fold degenerate levels, we obtain a total of 
2.704.156. Of this total, for seniority 5 = 0, i.e. no broken pairs, we have 924 states. Since the Hamiltonian does 
not connect states with different seniority 5, we can diagonalize a 924 x 924 matrix and obtain all eigenvalues with S 
= 0. Similarly, we can subdivide the Hamiltonian matrix into S = 2, S = 4,... and 5 = 12 (all pairs broken) blocks 
and obtain all 2.704.156 eigenvalues for a system with i = 12 levels and = 12 particles. As such, we have the 
exact density of levels and can compute observables like the entropy, heat capacity, etc. This numerically solvable 
model enables us to compute exactly the entropy in the microcanonical and the canonical ensembles for systems 
with odd and even numbers of particles. In addition, varying the relation 5 = d/G between the level spacing d and 
the pairing strength G may reveal features of the entropy that are similar to those of the experimentally extracted 
entropy discussed in the previous subsection. Recall that the experimental level densities represent both even-even 
and even-odd nucleon systems. 

Here we study two systems in order to extract differences between odd and even systems, namely by fixing the 
number of doubly degenerated single-particle levels to L = 12, whereas the number of particles is set to A'' = 11 and 
A^ = 12. 

These two systems result in a total of ~ 3 x 10^ eigenstates. In the calculations, we choose a single-particle level 
spacing of d = 0.1 MeV, which is close to what is expected for rare-earth nuclei. We select three values of the pairing 
strength, namely G = 1, 0.2, and 0.01, {5 = d/G = 0.1, 5 = d/G = 0.5, and 5 = d/G = 10), respectively. The first 
case represents a strong pairing case, with almost degenerate single-particle levels. The second is an intermediate case 
where the level spacing is of the order of the pairing strength, while the last case results in a weak pairing case. As 
shown below, the results for the latter resemble to a certain extent those for an ideal gas. 



1. Entropy 

The numerical procedure is rather straightforward. First we diagonalize the large Hamiltonian matrix (which is 
subdivided into seniority blocks) and obtain all eigenvalues E for the odd and even particle case. This defines also 
the density of levels p{E), the partition function, and the entropy in the microcanonical ensemble. Thereafter, we can 
obtain the partition function Z{T) in the canonical ensemble through Eq. ([76|). The partition function Z{T) enables 
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us, in turn, to compute the entropy S{T) by 

SiT) = kBlnZ{T) + {E{T))/T. (83) 
Since this is a model with a finite number of levels and particles, unless a certain smo othing is done, the microc anonical 



entropy may vary strongly from energy to energy (see for example the discussion in (Guttormsen et al, 2000)). Thus, 
rather than performing a certain smoothing of the results in the microcanonical ensemble, we will choose to present 
further results for the entropy in the canonical ensemble, using the Laplace transform of Eq. (^6|) . 

The results for the entropy in the canonical ensemble as functions of T for the above three sets of (5 = d/G are 
shown in Fig. For the two cases with strong pairing, we see a clear difference in entropy between the odd and the 
even system. The difference in entropy between the odd and even systems can be easily understood from the fact that 
the lowest-lying states in the odd system involve simply the excitation of one single-particle to the first unoccupied 
single-particle state and is interpreted as a single-quasiparticle state. These states are rather close in energy to the 
ground state and explain why the entropy for the odd system has a finite value already at low temperatures. Higher- 
lying excited states include also breaking of pairs and can be described as three-, five-, and more-quasiparticle states. 
For S — 10, the odd and even systems merge together already at low temperatures, indicating that pairing correlations 
play a negligible role. For a small single-particle spacing, also the difference in energy between the first excited state 
and the ground state for the odd system is rather small. 

For 5 — 0.5, we note that at a temperature of fc^T ^ 0.5 — 0.6, the even and odd systems approach each other^. The 
temperature where this occurs corresponds to an excitation energy (E) in the canonical ensemble of (E) ~ 4.7 — 5.0. 
This corresponds to excitation energies where we have 4 — 6 quasiparticles, seniority 5 = 4 — 6, in the even system 



and 5 — 7 quasiparticles, seniority 5 = 5 — 7, in the odd system (see for example the discussion in ( Guttormsen et al._ 
2000| )). For the two cases with strong pairing {S = 0.1 and S = 0.5), Fig. 30 tells us that at temperatures where we 



have 4 — 6 quasiparticles in the even system and 5 — 7 quasiparticles in the odd system, the odd and even system 
tend to merge together. This reflects the fact that pairing correlations tend to be less important as we approach 
the non-interacting case. In a simple model with just pairing interactions, it is thus easy to see where, at given 
temperatures and excitation energies, certain degrees of freedom prevail. For the experimental results, this may not 
be the case since the interaction between nucleons is much more complicated. The hope, however, is that pairing may 
dominate at low excitation energies and that the physics behind the features seen in Fig. ^ is qualitatively similar 
to the experimental information conveyed in Fig. Pq. 



2. The free energy 

We may also investigate the free energy of the system. Using the density of states, we can define the free energy 
F{E) in the microcanonical ensemble at a fixed temperature T (actually an expectation value in this ensemble), 

F{E) = -Tin [nNiE)e-f^^] . (84) 

Note that here we include only configurations at a particular E. 

The free energy was used by Lee and Kosterlitz ( Lee and Kostcrli"t^ , |1990 , 1991), based on the hi s togram approach 



for studying phase transitions developed by Ferrenberg and Swendsen ( pcrrenberg and Swendsen , |l9884 |b|) in their 



studies of phase transitions of classical spin systems. If a phase transition is present, a plot of F{E) versus E will 
show two local minima which correspond to configurations that are characteristic of the high and low temperature 
phases. At the transition temperature Tp, the value of F{E) at the two minima equal, while at temperatures below 
Tc, the low-energy minimum is the absolute minimum. At temperatures above Tc, the high-energy minimum is the 
largest. If there is no phase transition, the system develops only one minimum for all temperatures. Since we are 
dealing with finite systems, we can study the development of the two minima as a function of the size of the system 
and thereby extract information about the nature of the phase transition. If we are dealing with a second-order phase 
transition, the behavior of F{E) does not change dramatically as the size of the system increases. However, if the 
transition is first order, the difference in free energy, i.e., the distance between the maximum and minimum values, 
will increase with increasing size of the system. 

We calculate exactly the free energy F{E) of Eq. (|8^) through diagonalization of the pairing Hamiltonian of Eq. ( [2l| ) 
for systems with up to 16 particles in 16 doubly degenerate levels, yielding a total of ~ 4 x 10^ configurations. The 
density of states hence defines the microcanonical partition function. 



If we wish to make contact with experiment, we could assign units of MeV to fc^T and E. 
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For d/G = 0.5 and 16 single-particle levels, the calculations yield two clear minima for the free energy. This is 
seen in Fig. ^ where we show the free energy as a function of excitation energy using Eq. (|8^ ) at temperatures 
T — 0.5, T — 0.85, and T = 1.0. The first minimum corresponds to the case where we break one pair. The second 
and third minima correspond to cases where two and three pairs are broken, respectively. When two pairs are broken, 
corresponding to seniority 5 = 4, the free energy minimum is made up of contributions from states with 5 = 0, 2, 4. 
It is this contribution from states with lower seniority which contributes to the lowering of the free energy of the 
second minimum. Similarly, with three pairs broken, we have a free energy minimum which receives contributions 
from iS = 0, 2, 4, 6, yielding a new minimum. At higher excitation energies, population inversion takes place, and our 
model is no longer realistic. 

We note that for T = 0.5, the minima at lower excitation energies are favored. At T = 1.0, the higher energy 
phase (more broken pairs) is favored. We see also that at T = 0.85 for our system with 16 single-particle states 
and d/G = 0.5, the free-energy minima where we break two and three pairs equal. Where two minima coexist, we 
may have an indication of a phase transition. Note, however, that this is not a phase transition in the ordinary 
thermodynamical sense. There is no abrupt transition from a purely paired phase to a nonpaired phase. Instead, 
our system develops several such intermediate steps where different numbers of broken pairs can coexist. At e.g., 
T = 0.95, we find again two equal minima. For this case, seniority 5 = 6 and 5 = 8 yield two equal minima. This 
picture repeats itself for higher seniority and higher temperatures. 

If we then focus on the second and third minima, i.e., where we break two and three pairs, respectively, the difference 
Ai^ between the minimum and the maximum of the free energy can aid us in distinguishing between a first-order 
and a second-order phase transition. If AF/N, with N being the number of particles present, remains constant as 
N increases, we have a second-order transition. An increasing AF/N is, in turn, an indication of a first-order phase 
transition. It is worth noting that the features seen in Fig. |3^ apply to the cases with N = 10, 12, and 14 as well, with 
T = 0.85 being the temperature where the second and third minima equal. This means that the temperature where 
the transition is meant to take place remains stable as a functio n of number of single-pa r ticle levels and particles. 
This is in agreement with the simulations of Lee and Kosterlitz (Lee and Kosterlitz, 199C, 1991). We find a similar 
result for the minima developed at T = 0.95, where both 5 = 6 and 5 = 8. However, due to population inversion, 
these minima are only seen clearly for A'^ = 12, 14, and 16 particles. In Table II] we display AF/N for N = 10, 12, 
14, and 16 at T = 0.85 MeV. 

Table III reveals that AF/N is nearly constant, with AF/N « 0.5 MeV, indicating a transition of second order. 
This result is in agreement with what is expected for an infinite system. 

Before proceeding to the next method for classifying a phase transition in a finite system, we note the important 
result that for d/G > 1.5, our free energy, for N < 16, develops only one minimum for all temperatures. That is, 
for larger single-particle spacings, there is no sign of a phase transition. This means that there is a critical relation 
between d and G for the appearance of a phase transition- like behavior, being a reminiscence of the thermodynamical 
limit. This agrees also with e.g., the results for ultrasmall metallic grains (von Delft and Ralph, 2001). 



3. Distribution of zeros of the partition function 



Another way to classify the thermal behavior of finite systems requires the analytic continuation of the partition 



funct ion to the complex plane. Grossmann et aZ. ( Grossmann and Lchmann 



1969: Grossmann and Rosenhaucr, 1967 



1969) first introduced this technique for infinite systems. In these early works, the authors were able to indicate the 
nature of phase transitions by studying the densi ty of zeros (DOZ ) of the partition function. Borrmann et al. recently 
extended this idea to finite many-body systems (Borrmann et al, 2000| ). We implement the method by extending the 
inverse temperature to the complex plane (3 B = [3 + it. The partition function is then given by 



Z[B) = j dEp{E) exp{~^BE). 



B 



(85) 

Pk + iTk {k — I, - ■ ■ , A/") are complex 



Since the partition function is an integral function, the zeros Bk 
conjugated. 

Different phases are represented by regions of holomorphy that are separated by zeros of the partition function. 
These zeros typically lie on lines in the complex temperature plane. For a finite system, the zeros do not fall exactly 
on lines (they can be quite distinguishable depending on the size of the system) , and therefore the separation between 
two phases is more blurred than in an infinite system. The distribution of zeros contains the complete thermodynamic 
information about the system, and all thermodynamic properties are derivable from it. For example, in the complex 
plane, we define the specific heat as 



(86) 
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Hence, the zeros of the partition function become poles of Cv{B). A pole approaching the real axis infinitely closely 
causes a divergence at a real critical temperature Tc- The contribution of a zero Bk to the specific heat decreases 
with increasing imaginary part , so that thermodynamic properties of a system are governed by the zeros of Z lying 
close to the real axis. 

The distribution of zeros close to the real axis is approximately described by three parameters. Two of these 
parameters reflect the order of the phase transition, while the third indicates the size of the system. Let us assume 
that the zeros lie on a line. We label the zeros according to their closeness to the real axis. Thus n reflects the 
discreteness of the system. The density of zeros for a given Tk is given by 

^ ^^'^ ^ I {\B,\.,\ + \bJ-B,\) ' 

with A: = 2, 3, 4, • • •. A simple power law describes the density of zeros for small r, namely (/'(r) ~ r". If we use only 
the first three zeros, then a is given by 

In0(r3) - In0(r2) 

a ^ 



In T3 - In T2 



The final parameter that describes the distribution of zeros is given by 7 = tani^ ^ {P2 ~ Pi)I{t2 — ti). 

In the thermodynamic limit ti — > in which case the parameters a and 7 coincide with the infinite system limits 
discussed by Grossman et al. rossmann and Lchmann| , |1969| ; prossmann and Roscnhauci] , |1967| , |1969| ). For the 



infinite system, a = and 7 = yield a first-order phase transition, while for < a < 1 and 7 = or 7 ^ indicates 
a second-order transition. For arbitrary 7 third-order transitions occur when 1 < a < 2. For systems approaching 
infinite particle number, a cannot be smaller than zero since this causes a divergence of the internal energy. In small 
systems with finite ri, a < is possible. 

Continuation of the partition function to the complex plane is best interpreted by invoking a quantum-mechanical 
interpretation, namely 

Zip + ir) = fvA exp (^-irH^j cxp , (89) 

where the quantum-mechanical trace of an operator, projected on a specified particle number is given by 

TrAi^J2(a\pAi\a) , (90) 

a 

Pa is the number projection operator, and a runs over all many-body states. Since P represents the inverse tempera- 
ture, the thermally averaged many-body state is a linear combination of many-body states weighted by a Boltzmann 
factor, I vE'(/3,t = 0)) = exp(— /3i?Q,) | a), so that the partition function may be compactly written as 

ZiP + it) = i = 0) I vI/(/3, t = T)) . (91) 

Thus the zeros represent those times for which the overlap of the initial canonical state with the time-evolved state 
vanishes. In the r direction, the zeros represent a memory boundary for the system. 

In Fig. |3^ we show contour plots of the specific heat | Cy{B) \ in the complex temperature plane for iV = 11 (a), 14 
(b), and 16 (c) particles at normal pairing d/G = 0.5 and the = 14 (d) in the strong pairing limit, d/G = 1.5. The 
poles are at the center of the dark contour regions. We see evidence of two phases in these systems. The first phase, 
labeled / in Fig. |3^, is a mixed seniority phase, while the second phase, //, is a paired phase with zero seniority and 
exists only in even-A^ systems. No paired phase exists in the A^ = 11 system, and no clear boundaries are evident 
in the strong pairing case. We find that for (b) and (c) the DOZ are apparently distributed along two lines where 
the intersection occurs at ri, which is the point closest to the real axis. As the pairing branch (for (3 > (3i) only 
encompasses two points, we are unable to precisely determine a along this branch while 7 > 0. Based on our free 
energy results discussed above, we believe a along this branch will be positive. In the mixed phase branch (for (3 < (3i) 
we find 7 < 0, and a < in all normal-pairing cases. 



D. Level densities from shell-model Monte Carlo calculations 



We also a pplied shell-model Monte Carlo (SMMC) techniques to survey rare-earth nuclei in the Dy region ( White 
et a[], pOOO ). The goal of this extensive study was to examine how the phenomenologically motivated "pairing- plus- 
quadrupole" interaction compares in exact shell-model solutions with other methods. We also examined how the 
shell-model solutions compare with experimental data. 
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We discuss here one particular aspect of that work, namely level density calculations. Details may be found in 



( White et al , 2000 ). We used the Kumar-Baranger Hamiltonian with parameters appropriate for this region. Our 
single-particle space included the 50-82 subshell for the protons and the 82-126 shell for the neutrons. While several 
interesting aspects of these systems were studied in SMMC, we limit our discussion here to the level densities obtained 
for Dy. 

SMMC is an excellent way to calculate level densities. E{f3) is calculated for many values of /3 which determine the 
partition function, Z, as 

HZ{l3)/zm^- f d(3'E(P') . (92) 

Z{Q) is the total number of available states in the space. The level density is then computed as an inverse Laplace 
transform of Z . Here, the last step is performed with a saddle point approximation with (3~'^C = —dE/d/S: 

S{E) = f3E + lnZ{P) , (93) 
p(S) = (2^/3-2c)-i/2exp(5) . (94) 



The comparison of SMMC density in ^^^Dy with the Tveter et al. ( Tveter et al. , 1996 ) data is displayed in Fig. |3^. 
The experimental method can reveal fine structure, but does not determine the absolute density magnitude. The 
SMMC calculation is scaled by a factor to facilitate comparison. In this case, the factor has been chosen to make the 
curves agree at lower excitation energies. From 1-3 MeV, the agreement is very good. From 3-5 MeV, the SMMC 
density increases more rapidly than the data. This deviation from the data cannot be accounted for by statistical 
errors in either the calculation or measurement. Near 6 MeV, the measured density briefly flattens before increasing 
and this also appears in the calculation, but the measurement errors are larger at that point. 

The measured density includes all states from the theoretical calculation plus some others, so that one would expect 
the measured density to be greater than or equal to the calculated density and never smaller. We may have instead 
chosen our constant to match the densities for moderate excitations and let the measured density be higher than the 
SMMC density for lower energies (1-3 MeV). Comparing structure between SMMC and data is difficult for the lowest 
energies due to statistical errors in the calculation and comparison at the upper range of the SMMC calculation, i.e., 
E Ki lb MeV is unfortunately impossible since the data only extend to about 8 MeV excitation energy. 



VI. CONCLUSIONS AND OUTLOOK 

Pairing is an essential feature of nuclear systems, with several interesting and unsettled theoretical and experimental 
consequences, such as superfluidity and neutrino emission in neutron stars or pairing transitions in finite nuclei. 

This review is by no means an exhaustive overview; rather, our focus has been on the link between the nuclear 
many-body problem, and the underlying features of the nuclear force, and selected experimental interpretations and 
manifestations of pairing in nuclear systems. Our preferred many-body tools have been the nuclear shell model with 
its effective interactions and various many-body approaches to infinite matter. The common starting point for all 
these many-body approaches is, however, the free nucleon-nucleon interaction. 

Within this setting, we have tried to present and expose several features of pairing correlations in nuclear systems. 
In particular, we have shown that in neutron star matter (Sec. ||), pairing and superfluidity is synonymous with 
singlet ^Sq and triplet pairing up to densities 2-3 times nuclear matter saturation density. For singlet pairing 
it is the central part of the nucleon-nucleon force which matters, which within a meson-exchange picture can be 
described in terms of multi-pion exchanges. For triplet pairing, it is the two-body spin-orbit force which provides the 
attraction necessary for creating a positive pairing gap. Hyperon pairing, especially S" pairing, is also very likely. 
However, the actual size of these nucleon and/or hyperon pairing gaps in infinite neutron star matter is an unsettled 
problem and awaits further theoretical studies. A proper treatment of both short-range and long-range correlations is 
central to this problem. It will have significant consequences on the emissivity of neutrinos in a neutron star. Color 
superconductivity in the interior of such compact objects is also an entirely open topic. A similarly unsettled issue is 
the size of the triplet '^Si gap in symmetric matter or asymmetric nuclear matter. 



The above partial waves are also important for our understanding of pairing properties in finite nuclei. In Sec. HI 
we showed, for example, that the near-constancy of the excitation energy between the ground state with J = and 
the first excited state with J = 2 for the tin isotopes from ^^^Sn to ^'^^Sn is essentially due to the same partial waves 
which yield a finite pairing gap in neutron star matter. Moreover, a seniority analysis of the pairing content of the 
wave functions for these states shows that we can very well approximate the ground state with a seniority 5 = state 
(no broken pairs) and the first excited state in terms of a seniority 5 = 2 state (one broken pair). 
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Furthermore, we have used results from large-scale shell-model Monte Carlo and diagonalization calculations to 
extract information about isoscalar and isovector pairing and thermal response for /p-shell nuclei. One key result 
here was the decrease of T = 1 pairing correlations as a function of increasing temperature (up to about 1 MeV) 
and a commensurate buildup of structure in the specific heat curves at the sa me temperature. Information about 
proton-neutron pairing and the Wigncr energy has also been presented in Sec. 



Ill 



The important result here was 

that all J channels of the i nteraction contribute to the Wigne r energy, and that the J = 1 and J = Jmax channels 



contributed most (see also (Poves and Martincz-Pincdo 



elusive aspect of the nuclear pairing problem 
analysis. Our results from Sec. IV may indic a te th a t the T 



199^ )). Proton-neutron pairing is, however, a much more 
Its actual size, as is also the case in infinite matter, needs further 



role. For more information, see (Volya et al, 2002; [Zelevinsky et ai , 1996| ) 



part of the pairing interaction does not play a crucial 



Finally, we have tried to analyze recent experimental data on nuclear level densities in terms of pairing correlations. 
This was done in Sec. |^. These data reveal structures in the level density of rare-earth nuclei that can be interpreted 
as a gradual breaking of pairs. The experimental level densities can also be used to compute the thermal properties 
such as the entropy or the specific heat. The even systems exhibit a clear bump in the heat capacity. The temperature 
where this bump appears can be interpreted as a critical temperature for the quenching of pairing correlations. Similar 
features were also noted in Sec. Ill (see especially Figs. pO| and pl|) . More information was also obtained by studying 
the experimental entropy for even and odd nuclei with those extracted from a simple pairing model with a given 
number of particles and number of doubly-degenerate particle levels. We showed in Sec. ^ that the entropy of the 
odd and even system merge at a temperature which corresponds to the observed bumps in the heat capacity. This 
temperature occurs typically where we have 2-3 broken pairs. 

Within the framework of this simple pairing model, we showed also that for a finite system there is no sudden and 
abrupt transition to another phase, as we have in an infinite system. Rather, there is a gradual breaking of pairs as 
temperature increases. However, studying systems with different numbers of particles and levels, we presented also 
two possible methods for classifying the order of the transition. In order to perform these studies, we needed all 
eigenvalues from the simple pairing model in order to compute thermodynamical properties. Albeit there have been 
several interestin g theoretica l deve l opments of the solu t ion of the model Hamiltonian o f Eq. (^2|) or related models, 



see for example ( |Richardson| , 2002; Roman et al, ^002 ; Volya, 2002; Volya et al, 2001), we would like to stress that 



the investigation of thermodynamic properties requires a knowledge of all eigenvalues. 

Furthermore, an obvious deficiency of this simple model in nuclear physics is the lack of long-range correlations, 
which could, e.g., be expressed via quadrupole terms. A pairing-plus-quadrupole model, as discussed in Sec. Ill, would 
however, spoil the simple block-diagonalization feature in terms of seniority as a good quantum number. Such a model 
is however necessary, since the nuclear force is particular in the sense that the ranges of its short-range and long-range 
parts are rather similar. This means that short-range contributions arising from e.g., strongly paired particle-particle 
terms and long-range terms from particle-hole excitations are central for a correct many-body description of nuclear 
systems, from nuclear matter to finite nuclei. The difficulty connected with these aspects of the nuclear force means 
that further analysis of the thermodynamics of e.g., rare-earth nuclei can presently only be done in terms of large-scale 
shell-model Monte Carlo methods. 

It should be fairly obvious that we have only been able to cover a very limited aspect of pairing in nuclear systems. 
We have limited the attention to st a ble sy stems. However, pairing correlations are expected to play a special role in 

1996| ). There is currently a considerable experimental effort in nuclear physics. 



drip-line nuclei (Dobaczewski et al 



especially due to the advances from radioactive-beam and heavy-ion facilities, which have stimulated an exploration 
of nuclei far from stability. Many of these nuclei are weakly bound systems. Hence, due to strong surface effects, 
the properties of such nuclei are perfect laboratories for studies of the density dependence of pairing interactions. 
An experimental observable that may probe pairing corre l ation s is the pair transfer factor, which is directly related 
to the pairing density (see (Dobaczweski and Nazarewicz, 1998) for more details). The difference in the asymptotic 
behavior of the single-particle density and the pair density in a weakly bound system can be probed by comparing 
the energy dependence of one-particle and two-particle transfer cross sections. Such measurements, when performed 
on both stable and neutron rich nuclei, can hopefully shed some light on the asymptotic beha vior of pairing. An 
interesting system here is the chain of tin isotopes beyond ^'^^Sn. Various mean-field calculations (Dobaczewski et 



1996|) indicate that there i s a considerable increase in the pair transfer form factors for nuclei between ^^"Sn and ^''■^Sn 



TTTc 



1996| ). At the moment of writing, /3-decay properties of nuclei like ^^®Sn have just been studied 



( Dobaczewski et al, 
( [Shergur et al\ , |2002| ) 

From a many-body point of view, a correct treatment of these weakly bound systems entails an approach which 
encompasses a proper description of bound states and eventually features from the continuum. Such calculations have 
recent ly been mou nted within the framework of mean-field and Hartree-Fock-Bogoliubov (HFB) models (see (Grasso 
et a I, ^001| , |2002D ). A finite -range pairing interaction was included explicitely in the calculations. We mention here 
that the pairing terms in such me an-field calculations can be parameterized from microscopic many-body calculations, 
as demonstrated by Smerzi et al ( Smerzi et al , 1997). However, to include the continuum in a many-body description 
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such as the shell-model with an appropriate effective interaction is highly non-trivial. Even the determination of the 
effective two-body interaction is an open problem. The low-density studies of singlet ^Sq pairing in dilute Fermi 



systems (Heiselberg et al, 2000) clearly demonstrate that polarization terms cannot be neglected. 



We conclude by pointing to the strong similarities between pairing in the nuclear many-body p roblem and pairing in 



systems of trapped fcrmions (sec (Bruun and Mottclson, 2001; Heiselberg and Mottelson, 2002| ) for recent examples) 



There are also very strong cou plings to the experimental and theoretical developments of our understanding of ultra- 
small s uperconducting grains (Balian et al , 1999 ; Black et ai , 1996, 1997; von Delft and Ralphj , ^00 1| ; pukelsky and] 
SierQ [1999| ; [Mastellone and Falci] , |1998| ; [Ralph et al\ , |1995D . 
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TABLE I 2]*" — 0^ excitation energy for the even tin isotopes ^^^Sn for various approaches to the effective interaction. See 
text for further details. Energies are given in MeV. 

"«Sn "^Sn i^°Sn ^^^Sn ^^^Sn ^^^Sn ^^^Sn ""°S^ 

Expt 1.29 1.23 1.17 1.14 1.13 1.14 1.17 1.23 

Feff 1-17 1.15 1.14 1.15 1.14 1.21 1.28 1.46 

G-matrix 1.14 1.12 1.07 0.99 0.99 0.98 0.98 0.97 

^So G-matrix 1.38 1.36 1.34 1.30 1.25 1.21 1.19 1.18 

No ^So & ^P2 in G 0.15 -0.32 0.02 -0.21 



TABLE II Seniority v = overlap (first row) K^Sn; 0^|(5^) 2 [0)p and the seniority v = 2 overlaps (remaining rows) 

\{^Sn; Jf\D]jMiS'')^''\0)f for the lowest-lying eigenstates of '^^-^^°Sn. 
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TABLE III AF/N for T = 



0.85 MeV. See text for further details. 

N 10 12 14 16 

AF/N [MeV] 0.531 0.505 0.501 0.495 
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~ 10 km 



FIG. 1 Possible structure of a neutron star. 




FIG. 2 Plot of the NN interaction for the ^So channel employing the Argonne Vig interaction (Wiringa et ai, L995) 
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FIG. 3 Phase shifts for J < 2 partial waves as function of the incoming energy of two T = 1 nucleons in the lab system. 
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FIG. 4 P2 phase-shift p r edict ions of different potentials up to Eiah = 1-1 GeV, compared with the phase shift analysis of 
Arndt et al. (|Arndt et ai\, |l997|). 
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FIG. 5 Plot of the Argonne Vi8 ( Wiringa et al. , 1995) '^P2 partial wave contribution with and without the spin-orbit contribution. 
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FIG. 6 energy gap in neutron matter with the CD-Bonn, Nijmegen I and Nijmegen II potentials. In addition, we s how the 
results obtained fro m ph ase shifts only, Eqs. (^)-(^3|), and the effective range approximation of Eq. (jsH]). Taken from (Elgar0y 
and Hjorth-Jenser , 199^). 
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FIG. 7 Top panel: The angle-averaged ^P2-^F2 gap in neutron matter depending on the Fermi momentum, evaluated with free 
single-particle spectrum and different nucleon-nucleon potentials. Middle panel: The gap evalu ated with BHF spec tra. Bottom 
panel: The gap with the CD-Bonn potential in different approximation schemes. Taken from (Baldo et al, 1998b). 
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Figure 8: dean_hjorthjensen_fig8.ps 



FIG. 8 Effective masses derived from various interactions in the BHF approach. Taken from (Baldo et at, 19981: ) 
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FIG. 9 Upper box: Proton pairing in / 3-stable matte r for th e ^5*0 partial wave. Lower box: Neutron pairing in /3-stable matter 
for the partial wave. Taken from (Elgar0y et al., 1996b). 
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Figure 10: dean_hjorthjensen_figl0.ps 



FIG. 10 EOS for sym metric nuclear matter with the NN potentials and many-body methods described in the text. Taken from 
(|Elgar0y et~^, |l998|). 
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Figure 11: dean_hjorthjensen_figll.ps 



FIG. 11 ^Si-lDj_ energy gap in nu clear matter, calculated in non-relativistic (full lines) and relativistic (dashed line) approaches. 
Taken from (|Elgar0y et a;.i|l99g|). 
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FIG. 12 Energy gap in different approximations for tfie self-energy. Tfie upper line (dashes with crosses) stands for the free 
single-particle spectrum with a standard BCS approach while the upper solid line arises from the BHF approach of Eq. (^) 
and the standard BCS approach. The lower lines arise from solving Eq. (^) for the pa i ring gap with differe nt approaches to 
the se l f-energy, for further details se e (Lombardo et al, 2001b; Lombardo and Schulze, 2001). Taken from ( [Lombardo et al 
2001b|; [Lombardo and Schulz4 |200l|). 
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FIG. 13 Diagram (a) is the second-order diagram with particle-hole intermediate states. The external legs can be particles or 
holes. Diagrams (b) and (c) are examples of third-order TDA or RPA diagrams. The dotted lines represent the interaction 
vertex. 
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FIG. 14 The ^Sn gap in pure neut ron matter predicted in several publications taking account of polarization effects. Taken 
from ( Lombardo and Scliulz4 2001) 
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FIG. 15 SMMC expect ation values ol proton and neutron BCS-like pairs after subtraction of the Fermi gas value. Taken from 
(Langanke et al., 1995). 
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FIG. 16 Proton and neutron occupation numbers for various chains of isotopes. Taken from (Langanke et al., 1995) 
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FIG. 17 Content of isovector 0"*" pairs and isoscalar 1^ pairs in the ground states of the isotopes ^*Fe. The upper panel 
shows the diagonal matrix elements of the pair matrix Mac- The index a = 1, ..,4 refers to 0^ pairs in the /7/2, P3/2, f5/2, 
and pi/2 orbitals, respectively. For the isoscalar pairs a = 1,2,3 refers to {fT/2)^, (/7/2/5/2), and pairs, respectively. 

The middle panel gives the eigenvalues of the pair matrix; for the isoscalar pairs, only the 3 largest are s hown. The lower pane l 
gives the eigenvalues of the pair matrix for the uncorrelated Fermi gas case using Eq. (|63L Taken from (Langanke et ai, 1996). 



J-O, T-1 pn 




46 48 50 52 54 56 58 60 



50 52 54 56 58 60 



FIG. 18 Largest eigenva lue for the various isovector 0"'" pairs in the N — Z nuclei in the mass region A — 4S — 56. Taken from 
(Langanke et al, 1997a). 
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FIG. 19 Calculated displacement ew of the binding energy of 2^Mg and ^*Cr from the average parabolic (A'^ — Z)'^ behavior 
along an isobaric chain. Shell-model calculations were performed in the QTiuj configuration space. The results of calculations 
for the binding energies of even-even nuclei along the A=48 chain (normalized to ''^Cr) are shown in the insert. The values of 
ew were obtained using the shell-model Hamiltonian with the J = 1, 2, .., Jmax, T = matrix elements removed. For instance, 
the result for Jniax=3 corresponds to the variant of calculations in which all the two-body matrix elements between states 
|jij2>/r=0) with J=l,2,3, were put to zero. The results are normalized to th e full shell-rnodel v alue e^*''' (Jmax=0). The 



Coulomb contribution to the binding energy has been disregarded. Taken from (Satula et al.. 1997) 
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FIG. 20 Temperature dependence of various observables in ^''Fe. Monte Carlo points with statistical errors are shown at 
each temperature T. In the left-hand column, the internal energy, U, is calculated as (H) — Eo, where H is the many-body 
Hamiltonian and Eo the ground state energy. The heat capacity C is calculated by a finite-difference approximation to dU/dT, 
after U{T) has been subjected to a three-point smoothing, and the level density parameter is a = C/2T. In the right-hand 
column, we show the expectation values of the squares of the proton and neutron BCS pairing fields. For comparison, the 
pairing fields calculated in an uncorrelated Fermi gas are shown by the solid curve. The moment of inertia is obta ined from 
the expectation values of the square of the total angular momentum by 7 = f3{j^)/3. Taken from (Dean et al, 1994). 
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FIG. 21 The heat capacity of even-even (left panel) and odd-even (right panel) iron isotopes. Taken from ( Liu and Alhassi'^ , 
200l|). 
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FIG. 22 Temperature dependence of various observables in ^°Mn. The left panels show (from top to bottom) the total energy, 
the moment of inertia, and the proton-neutron BCS pairing fields, while the right panels exhibit the expectation values of the 
isospin operator (T^), the isovector J = proton-proton BCS pairing fields, and the isoscalar J — 1 pairing strength, as defined 
in the text. For comparison, the solid lines indicate the Fermi gas values of the BCS pairing fields and J = 1 pairing strength. 
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FIG. 23 Pair correlations, as defined in Eg. pq, for isovector 0^ and isoscalar 1^ pairs for ^^Fe and ^®Cr, as functions of 



temperature. Taken from (Langanke et ai, 1996) 



FIG. 25 Distribution of the pair transfer cocfRciont <P,)|,>, pairing strength <k>, and deformation Q2 for the random 
interactions leading to a J'^ = 0"*" ground state in the shell-model description. From left to right we report the results for 
N — Z = 4, 2, and 0. The arrows indicate the number of results in a bin when it is out of scale. 
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FIG. 26 Observed entropy for ^'^^•^^^Dy as function of excitation energy E. Taken from ( Guttormsen et al. 



200C). 
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FIG. 27 Observed entropy for I'^i-i^^yj-, g^g function of excitation energy E. Taken from (Guttormsen et al. 



200C). 
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FIG. 28 Semi-experimental entropy S for 
perature fcsT. Taken from (Guttormsen et al 



pDv and 

iooq). 



Yb calculated in the canonical ensemble as a function of tem- 
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FIG. 29 Semi-experimental heat capacity as a function of temperature (left panels) and energy {E) (right panels) in the 
canonical ensemble for ^'^^''^^^Dy and ^^^'^^^Yb. The dashed lines describe the approximate Fermi gas heat capacity. The 
arrows indicate the first local maxima of the experimental curve relative to the Fermi gas estimates. The dashe d-dotted lines 
descr ibe extrapolated estimates of the critical temperature Tc. Tc is indicated by the vertical lines. Taken from ( ^chiller et al. , 
2OO1I ). 
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FIG. 30 Entropy in the canonical ensemble as a function of temperature ksT for odd and even systems for 5 = 0.1 (upper 
panel), 5 — 0.5 (central panel), and 5 = 10 (lower panel). If we wish to make contact wit h experiment, one can assign units of 
MeV to fcflT. The entropy S/ks is dimensionless. Taken from (Guttormsen et ai, 2000). 
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FIG. 31 Free energy from Eq. (||) at T = 0.5, 0.85 and T = 1.0 MeV with d/G = 0.5 with 16 particles in 16 doubly degenerate 
levels. All energies are in units of MeV and an energy bin of 10~^ MeV has been chosen. 
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FIG. 32 Contour plots of the specific heat in the complex temperature plane for a) A'^ = 11, b) A'^ = 14, and c) A*' = 16 particles. 
Panel d) shows the TV = 14 case with weak pairing. The spots indicate the locations of the zeros of the canonical partition 
function. 
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FIG. 33 SMMC density vs. experimental data in ^'^^Dy. 



